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Abstract 


It  has  long  been  assumed  that  the  normalized  spectrum  of  a  radiated  field 
remains  invariant  on  propagation.  Recent  developments  in  coherence  theory  have 
demonstrated  that,  in  general,  the  normalized  spectrum  of  a  radiated  field  changes  on 
propagation  depending  on  the  state  of  coherence  and  the  original  spectrum  of  the 
source.  In  this  thesis  we  examined  the  effect  of  the  spatial  correlations  of  the  source 
fluctuations  on  the  radiated  fields. 

We  first  considered  the  effect  of  source  correlations  on  the  total  radiated 
power  using  the  concept  of  "radiation  efficiency".  For  a  class  of  planar  sources, 
known  as  Schell-model  sources  we  showed  that  the  radiation  efficiency  and  hence  also 
the  total  emitted  power  increase  with  increasing  source  size  or  increasing  source 
correlation  length.  Interestingly  enough,  our  results  indicate  that  any  source  of  this 
class  whose  linear  dimensions  are  larger  than  about  a  wavelength  has  a  radiation 
efficiency  that  exceeds  90%. 

In  investigating  the  radiation  efficiency  of  three-dimensional,  partially 
coherent,  primary  sources  we  develop  a  method  for  finding  the  correlation  function  that 
maximizes  the  radiation  efficiency  and  the  total  power  emitted  by  such  sources.  In  the 
case  of  quasi-homogeneous  sources  we  show  that  the  optimal  degree  of  spatial 
coherence  is  smkr'Ikr’  where  r'  is  the  spatial  offset  and  k  is  the  wave  number.  The 
significance  of  this  result  is  discussed  in  connection  with  blackbody  radiation. 

In  considering  the  effects  of  the  source  correlation  on  the  spectrum  of  the 
radiation  we  analyze  a  simple  physical  configuration  of  two  small  sources.  Our 
calculations  demonstrate  that  one  can  choose  correlation  functions  that  give  rise  to  line 
narrowing,  line  broadening,  line  shifting  and  line  splitting.  Similar  results  are  also 
obtained  when  the  full  electromagnetic  nature  of  the  sources  is  taken  into  account.  We 
illustrate  this  fact  by  considering  radiation  from  two  partially  correlated  linear  dipoles. 


To  gain  a  different  perspective  for  the  effects  of  the  spatial  correlations  on  the 
spectrum  we  consider  an  example  in  coherent-mode  representation.  In  this  example  we 
introduce  the  concept  of  "spectral  modifier"  to  show  the  small  spectral  changes  that 
occur  on  propagation  of  single  coherent  modes,  and  we  then  examine  the  spectral 
effects  when  several  modes  are  present. 

In  the  final  part  of  the  thesis  we  consider  the  spectral  effects  arising  ft'om 
sowce  correlations  when  the  fields  propagate  in  homogeneous  or  inhomogeneous 
media.  Our  calculations  elucidate  the  development  of  spectral  shifts  as  a  function  of  the 
propagation  distance  and  illustrate  the  various  spectral  effects  for  light  propagation  in 
graded-index  fibers. 
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is  10°  and  the  source  parameters  are  JtoOl  =  20  and  kQOg  =  10.  Page  158. 

Figure  6.9:  Frequency  shift  Av  versus  the  propagation  distance  kQz  in  a 
dispersive  graded-index  medium  (a).  Curve  (b)  shows  Av  when  the  inhomogeneous 
nature  of  the  medium  is  ignored  by  setting  a  =  0.  Curve  (c)  shows  the  corresponding 
result  for  free-space  propagation.  The  observation  angle  is  10°  and  the  source 
parameters  were  chosen  to  be  =  20  and  itoOg  =  10.  Page  160. 

Figure  6.10:  Frequency  shift  Av  as  a  function  of  propagation  distance  in  a  graded- 
index  fiber  (solid  line).  The  frequency  shifts  arc  calculated  for  observation  at  a  fixed 
distance  10/^  from  the  center  of  the  fiber  and  =  20  and  kQC^  =  10.  The  dashed 
line  shows  the  frequency  shifts  when  the  frequency  dependence  of  a  is  ignored  by 
setting  [a((0o)/^o  =  0.00048].  Page  162. 
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Introduction 

The  term  spectrum  plays  an  important  role  in  most  areas  of  science.  In  optics  it  refers 
to  the  distribution  of  energy  as  a  function  of  frequency.  This  concept  has  deep  roots 
in  the  development  of  optical  physics.  The  first  reported  experiments  demonstrating 
that  sunlight  consists  of  multi-color  contributions  were  performed  by  Sir  Isaac 
Newton  and  published  as  early  as  1666.^  Since  the  days  of  Newton,  observations 
and  measurements  of  spectra  have  became  a  leading  research  tool  that  was  largely 
responsible  for  the  development  of  quantum  mechanics  and  has  made  significant 
contributions  to  areas  such  as  thermodynamics,  chemistry,  astronomy  and  metrology. 

In  spite  of  the  fact  that  relatively  accurate  measurements  of  spectra  were 
made  as  early  as  18172,  only  rough  guesses,  intuition  and  an  ample  measure  of  luck 
made  the  experimental  results  useful  for  the  development  of  the  various  sciences  that 
were  based  on  spectroscopy.^  Implicit  in  all  spectroscopic  measurements  is  the 
assumption  that  the  spectrum  of  the  field  measured  by  various  means  is  equal  to  the 
spectrum  of  the  source  generating  it,  even  when  the  light  has  propagated  a  significant 
distance  from  the  source  This  assumption,  as  was  demonstrated  recently,  is  valid 
only  in  special  cases.  In  general,  the  spectrum  of  partially  coherent  light  changes  on 
propagation.  The  extent  and  type  of  the  changes  are  determined  by  the  state  of 
coherence  of  the  source  and  separately  by  the  propagation  medium. 

Although  the  concept  of  coherence  was  recognized  many  years  ago,  it 
appears  that  scientists  have  not  considered  the  effects  of  coherence  on  the  spectrum 
until  recently.^  Experimentally,  it  is  known  that  the  most  common  sources  of  light  do 
indeed  produce  a  spectrum  that  does  not  vary  appreciably  on  propagation.  Although 
several  cases  of  spectral  changes  were  encountered,  they  were  either  attributed  to  the 
Doppler  effect^  or  left  as  an  experimental  uncertainty.^  Theoretically,  the  main 
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reason  that  coherence  effects  on  the  spectrum  have  not  been  investigated  earlier  may 
be  due  to  the  fact  that  for  many  years  the  main  quantities  used  in  coherence  theory 
were  space-time  correlations.  In  particular,  the  second  order  correlation  function 
most  widely  used  was  the  mutual  coherence  function.'^  One  of  the  major 
developments  of  statistical  optics  in  the  last  decade  was  the  formulation  of  coherence 
theory  in  the  space-frequency  domain,  which  shifted  the  emphasis  from  the  mutual 
coherence  function  to  the  cross-spectral  density  function.  The  transition  to  space- 
frequency  domain  simplifies  some  of  the  theoretical  techniques  used  to  solve 
boundary  value  problems  and  focuses  the  emphasis  on  the  frequency-dependence  of 
all  the  relevant  physical  quantities. 

Wolf  first  derived  conditions  under  which  the  spectrum  of  light  generated 
by  a  quasi-homogeneous  source  remains  invariant  after  propagation.^*’  According  to 
that  analysis,  planar  secondary  quasi-homogeneous  sources  whose  degree  of 
correlation  obeys  a  certain  scaling  law  produce  light  whose  spectrum  is  independent 
of  the  direction  of  observation  and  is  equal  to  the  source  spectrum.  Additional 
investigations  have  since  shown  the  effects  of  spatial  correlations  in  various  physical 
situations. 

The  simplest  physical  system  that  was  considered  was  that  of  two  small 
radiating  sources.®  The  simplicity  lies  in  the  fact  that  both  sources  are  assumed  to  be 
essentially  point  sources  and  the  spatial  correlation  is  a  function  that  depends  on  the 
position  of  the  two  sources  and  on  the  frequency.  It  has  been  shown  that  even  this 
simple  system  can  produce  spectral  line  shifts  as  well  as  actual  modulation  of  the 
radiated  spectrum.  These  theoretical  predictions  have  since  been  verified 
experimentally  by  several  groups.^ 

In  considering  different  spatial  correlations,  the  next  level  of  complexity  is 
that  of  extended  secondary  planar  sources.  Here  the  spatial  correlations  are  taken 
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between  any  two  points  in  the  source  plane.  Most  of  the  investigations  relating  to 
planar  sources  examined  the  spectral  shifts  occurring  on  the  propagation  of  light  from 
the  source  plane,  lo  Additional  woric  explored  more  general  spectral  changes  such  as 
line  broadening,  line  narrowing  and  line  distortion  in  partially  coherent  sources  whose 
cross-spectral  density  consists  of  several  Hermite-Gaussian  modes. A  number  of 
experiments  in  this  area  have  also  been  reported. 

When  the  radiating  sources  are  three-dimensional,  the  spatial  correlations 
have  additional  degrees  of  freedom  compared  with  the  spatial  correlations  of  planar 
sources. The  effects  of  the  spatial  correlations  on  the  spectrum  of  the  light 
produced  by  three-dimensional  sources  are,  however,  similar  to  those  produced  by 
planar  sources, 

The  investigation  of  changes  in  the  spectrum  of  partially  coherent  light  on 
propagation  have  dealt  so  far  mainly  with  propagation  in  free  space.  In  many 
practical  applications  light  propagates  in  dispersive  media  that  are  homogeneous  or 
inhomogeneous.  In  these  cases  it  is  important  to  estimate,  for  example,  the  changes 
in  the  spectrum  after  propagating  a  certain  distance  in  the  medium.  15 

The  spectrum  of  the  field  is  not  the  only  physical  quantity  affected  by  source 
correlations.  Qosely  related  quantities  are  the  radiated  power  at  a  given  frequency, 
the  directivity  of  a  radiating  system  and  the  shape  of  intensity  profile  of  light  beams 
produced  by  partially  coherent  sources.  The  effects  of  source  correlations  on  the 
spatial  distribution  of  the  radiated  intensity  have  been  investigated  for  several  years,  i^ 
Of  interest  in  this  thesis  are  the  radiation  efficiency,  the  directivity  and  the  shape  of  the 
intensity  profile  of  partially  coherent  optical  beams. 

The  radiation  efficiency  of  partially  coherent  sources  is  a  concept  that  has 
been  developed  as  a  measure  for  the  amount  of  light  generated  by  a  source  of  a  given 
intensity  profile  and  with  varying  forms  of  spatial  correlations.  Most  of  the  early 
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work  on  radiation  efficiency  involved  quasi-homogeneous  sources.^^  Other  types  of 
partially  coherent  sources  including  planar  and  three-dimensional  sources  have  since 
been  investigated.*"^  Considering  the  effects  of  the  correlation  ot  the  radiation 
efficiency,  it  is  natural  to  look  for  the  spatial  correlation  that  maxinoizes  the  radiation 
efficiency  and  hence  the  total  power  emitted  by  a  source  of  a  prescribed  intensity 
profile.** 

The  total  emitted  power  does  not  give  much  information  about  the 
distribution  of  the  radiation  in  space.  The  radiant  intensity  at  a  given  frequency  is 
important  for  the  characterization  of  the  directivity  of  fields  produced  by  partially 
coherent  sources  and  in  calculations  of  the  intensity  profiles  of  partially  coherent 
beams.  Since  the'radiant  intensity  depends  (mi  spatial  correlation,  so  do  the  directivity 
and  the  intensity  profile.  Of  particular  interest  in  this  respect  arc  the  intensity  profiles 
of  the  so-called  Gaussian  Schell-model  sources*^  that  are  closely  related  to  certain 
laser  radiators. 
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1.2  Organization  of  the  thesis 

This  thesis  describes  several  investigations  on  the  effect  of  source  correlation  on  the 
physical  properties  of  the  radiated  field.  Specifically,  we  discuss  the  dependence  of 
the  spectrum  of  light  produced  by  partially  coherent  sources  on  the  spatial  correlations 
of  the  source  fluctuations.  Most  of  the  material  included  in  the  following  chapters  has 
already  been  reported  in  papers  that  are  already  published  or  have  been  submitted  to 
publication.  The  pertinent  references  are  listed  in  Appendix  E.  Each  chapter  is 
concerned  with  a  specific  aspect  of  the  theory,  and  the  presentation  does  not  always 
follow  the  chronological  development  of  the  work.  Throughout  this  thesis  attempts 
have  been  made  to  thoroughly  document  the  relevant  literature.  End  notes  to  every 
chapter  contain  extensive  lists  of  references  in  addition  to  supplemental  statements  that 
could  not  be  incorporated  in  the  main  text  Since  much  of  the  research  concentrates 
on  closely  related  areas,  some  of  the  fundamental  papers  may  be  referenced  more  than 
once,  where  appropriate. 

In  the  second  part  of  this  chapter  we  review  some  elements  of  coherence 
theory  and  establish  the  notation  and  conventions  used  throughout  the  rest  of  the 
thesis.  Our  treatment  is  based  on  classical  statistical  optics  and  on  the  consideration 
of  statistically  stationary  fields.  Most  of  the  formal  treatment  is  based  on  space- 
frequency  representation  so  that  quantities  such  as  the  radiant  intensity,  the  total 
power  and  the  directivity  refer  to  their  respective  values  at  a  single  frequency.  We 
also  present  some  of  the  basic  relations  in  the  space-time  domain  in  order  to  illustrate 
the  advantages  of  using  the  space-frequency  representation.  Following  the  discussion 
of  the  cross-spectral  density  function  we  describe  model  sources  which  are  frequently 
used  in  the  theory  as  well  as  in  the  formulation  of  coherent-mode  representation. 


Introduction 


Chapua- 1 


7 


Chapter  2  contains  a  description  of  investigations  of  the  radiation  efficiency 
and  directivity  of  fields  produced  by  paniaiiy  coherent  planar  and  three-dimensional 
sources.  Specifically  we  examine  the  radiadon  efficiency  of  planar  Gaussian  Schell- 
model  sources.  We  determine  the  spadal  correladon  which  maximizes  the  radiadon 
efficiency  for  three-dimensional  sources  of  prescribed  intensity  profiles.  While  we 
show  that  the  maximum  radiadon  efficiency  of  planar  sources  is  obtained  when  the 
source  is  fully  coherent,  three-dimensional  coherent  sources  arc,  in  general,  highly 
inefficient  Since  the  radiadon  efficiency  and  the  direcdvity  arc  evaluated  at  a  single 
frequency,  it  is  necessary  to  account  for  the  original  source  spectrum  to  obtain  the 
total  efficiency  or  directivity.  Some  of  the  spectral  effects  are  considered  in  the 
following  chapters. 

In  chapter  3  we  begin  our  discussion  of  correladon  induced  spectral 
changes.  We  consider  a  simple  physical  configuradon  of  two  small  sources  and 
describe  some  of  the  possible  ways  in  which  the  spectrum  of  light  can  be  modulated 
by  appropriately  modifying  the  spadal  correladons.  The  analysis  up  to  this  point  is 
based  on  scalar  theory.  In  chapter  4  the  analysis  is  extended  to  the  full 
electromagnedc  fields,  by  considering  the  radiadon  produced  by  two  linear  electric 
dipoles.  This  chapter  expands  the  analysis  of  the  spectral  moduladon  technique  to 
include  direcdonal  effects  and  the  vector  properties  of  the  field.  A  comparison 
between  the  results  of  chapter  3  and  chapter  4  shows,  as  one  might  perh^s  expect, 
that  the  electromagnedc  treatment  reduces  to  the  results  of  scalar  theory  in  certain 
limits. 

The  theoretical  predictions  of  spectral  changes  in  chapters  3  and  4  are 
discussed  without  suggesting  any  possible  mechanism  that  can  produce  the  required 
spatial  correlations.  In  chapter  5  we  consider  the  radiation  from  a  secondary  planar 
source  whose  cross-spectral  density  consists  of  several  coherent  modes.  The  choice 
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of  the  Heimite-Gaussians  as  the  functional  form  of  the  coherent  inodes  makes  the 
analysis  relevant  to  applications  where  independent  laser  modes  are  present  in  a 
source.  The  theory  develojjcd  in  this  chapter  shows  the  spectral  changes  that  occur 
on  propagation  of  individual  coherent  modes.  When  the  radiation  consists  of  several 
coherent  modes  the  changes  in  the  spectrum  can  be  attributed  to  correlation  effects 
which  are  manifested  by  the  relative  strengths  of  the  individual  modes  present  in  the 
source. 

The  analysis  presented  in  chapters  2  through  5  concerned  propagation  of 
light  in  free  space.  In  most  practical  applicaticMis  light  passes  through  optical  elements 
and  through  media  whose  response  is  frequency-dependent  In  chapter  6  we  consider 
propagation  of  partially  coherent  light  through  homogeneous  and  inhomogeneous 
media,  as  well  as  the  effects  of  dispersion.  In  particular,  we  consider  propagation 
through  optical  fibers  whose  index  of  refraction  has  a  quadratic  profile  (so-called 
Selfoc  fibers).  As  one  may  expect,  spectral  changes  that  take  place  on  propagation  in 
free  space  are  somewhat  enhanced  by  propagation  in  a  medium  characterized  by  a 
frequency  dependent  index  of  refraction. 

Throughout  this  thesis  various  symbols  were  used  to  denote  numerous 
physical  and  mathematical  entities.  We  made  every  effort  to  define  all  symbols  the 
first  time  that  they  occur  in  the  text  In  Appendix  F  we  present  a  list  of  the  main 
symbols  and  their  definitions. 
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1.3  Elements  of  coherence  theory 

In  this  section  wc  summarize  some  basic  results  of  coherence  theory  which  are 
essential  to  the  development  of  the  material  presented  in  the  rest  of  the  thesis.  We 
focus  our  attention  on  the  development  of  space-frequency  representation,  and  we 
discuss  properties  of  the  cross-spectral  density  and  the  so-called  coherent-mode 
representation.  This  review  also  includes  short  derivations  of  the  mam  relations 
involving  various  physical  quantities  of  interest. 


The  concept  of  spatial  correlation 

Consider  an  aperture  in  an  opaque  screen^o  that  is  illuminated  by  a  light  source  as 
shown  in  Fig.  1.1.  The  field  distribution  in  the  plane  of  the  aperture  is  known  as  a 
secondary  source.  The  concept  of  a  secondary  source  is  useful  in  describing  practical 
situations  where  one  considers  the  light  incident  on  the  entrance  or  exit  pupil  of  an 
optical  system,  or  when  describing  a  particular  distribution  of  light  without  regard  to 
the  physical  system  which  produces  the  radiation  in  the  first  place.  We  denote  the 
field  strength  at  two  typical  points  in  the  plane  of  the  aperture  and  at  two  different 
times  by  complex  analytic  signals^i  y(Pj,  rj),  V'lFj.  respectively.  The 
correlation  between  the  fluctuations  of  the  field  at  the  two  points  is  characterized  by 
the  average 


C(/»,F2;r,,r2)  =  (v*(P,,  i^))  ■ 


(1.1) 
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Figure  1.1:  A  Secondary  source  consisting  of  an  aperture  in  a  planar  t^que 
screen.  and  P2  are  two  typical  points  in  the  plane  of  the  apmure. 

This  type  of  average,  which  occurs  frequently  in  our  analysis,  is  taken  over  an 

ensemble  of  realizations  of  similar  systems,  and  it  is  denoted  by  the  angular 

brackets.22 

In  this  document  we  consider  only  fields  which  arc  stationary  in  the  wide 
sense,  i.e.  fields  wliose  average  is  a  constant  which  is  independent  of  time,  and 
whose  two-point  correlation  F  depends  on  the  two  time  arguments  only  through  the 
difference  ^2  *  namely 

r(ri,  rj;  T)  =  (v*(r,,  t)Vir^,  t  +  x))  .  (1.2) 

The  space-time  correlation  I\ri,  ri,  t)  is  called  the  mutual  coherence  function?^ 
which  had  been  widely  used  together  with  its  normalized  form,  the  complex  degree  of 
coherence 
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y(»i,r2;T) 


rCr^.r^-.T) 


whose  magnitude  satisfies  the  constraint 

|y(r,.r2;T)|$l  . 


(1.3) 


(1.4) 


In  Eq.  (1.4)  the  lower  limit,  zero,  is  obtained  foi  a  temporally  incoherent  field  and  the 
upper  limit  ccnresponds  to  a  tempcrally  coherent  field. 

The  mutual  coherence  function  conveniently  facilitates  calculations 
involving  interference  of  partially  coherent  Ught.  For  calculations  involving  the 
propagation  of  partially  coherent  light  it  is  usually  better  to  use  the  cross-spectral 
density  function.  Before  we  define  it,  we  consider  the  generalized^^  Fourier 
representation  for  the  field  V(r,  t)  in  the  form 


V(r,r)  =  J v(r;<o)e-*“  da  .  (1.5) 

0 

The  cross-spectral  density  function  is  defined  by  the  ensemble  average 

Wiryr2;a)S{a-a')  =  (y{r^;a)v{r2;a')^  ,  (1.6) 

where  5((0  -  co’)  is  the  Dirac  delta  function.  Its  appearance  in  Eq.  (1.6)  arises  from 
the  assumption  of  stationarity  of  the  field.  I  is  important  to  observe  that  the  cross- 
spectral  density  function  defined  by  Eq.  (1.6)  is  a  measure  of  the  spatial  correlation  of 
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the  field  fluctuations  at  a  single  frequency  to.  Because  of  our  assumption  of 
stationarity,  the  cross-spectral  density  function  has  no  information  about  the  spatial 
correlation  of  the  fluctuations  of  the  field  at  different  temporal  frequencies.  As  we 
show  below,  the  fact  that  the  cross-spectral  density  is  a  spatial  correlation  at  a  single 
frequency  makes  it  particularly  suitable  for  calculations  where  the  frequency 
dependence  of  the  variables  is  important. 

The  cross-spectral  density  function  and  the  mutual  coherence  function  are 
related  by  the  generalization  of  the  well  known  Wiener-Khintchine  theorem^ 


W(ri,r2;n))=  dx  ,  (1.7) 

— oo 

ao 

r(r^,r^,t)  =  jw(r^,r2,co)e-‘^dv .  (1.8) 

0 

This  relatively  simple  relation  is  justified  on  the  basis  of  a  substantial  amount  of 
mathematical  analysis  that  is  outside  the  scope  of  this  manuscripL^^ 

Before  leaving  the  subject  of  the  mutual  ccherence  function  we  compare  the 
expressions  governing  the  propagation  of  the  mutual  coherence  function  and  of  the 
cross-spectral  density  function  to  the  far  zone.  This  comparison  illustrates  one  of  the 
advantages  of  using  the  cross-spectral  density  in  propagation  calculations.  In  free 
space,  the  mutual  coherence  function  is  given  by 


Intrcxluction 


Chapter  1 


Introduction 


Chapter  1 


14 


where  u-^  are  projections  of  the  unit  vectors  iij  on  the  source  plane^*,  and  is 

the  four-dimensional  Fourier  transform  of  the  cross-spectral  density  in  the  source 
plane 


(1.13) 


In  this  notation  k  in  Eq.  (1.12)  is  the  wavenumber  associated  with  frequency  co. 


k 


In 

T 


2ttv  _  (o 
c  c 


(1.14) 


and  X  is  the  wavelength  of  the  field. 

It  is  clear  from  Eqs.  (1.9)-(1.12)  that,  in  general,  it  is  simpler  to  use  the 
cross-spectral  density  function  when  calculating  the  propagation  of  partially  coherent 
fields  in  free  space.  In  the  rest  of  this  section  we  review  the  properties  of  the  cross- 
spectral  density  and  summarize  some  of  the  relationships  between  the  cross-spectral 
density  and  various  physically  measurable  quantities.^ 

The  first  measurable  quantity  that  we  consider  is  the  spectrum  5(r,  (o)  which 
is  simply  the  ‘diagonal’  component  of  the  cross-spectral  density,  i.e., 

5(r,<u)  =  H^(r,r;tt))  .  (1.15) 

Sir,  (o)  may  represent  the  spectrum  of  the  field  or  the  spectrum  of  the  source.  It  is 
usually  supplemented  by  a  subscript  or  superscript  giving  it  the  proper  designation, 
e.g.,  5(“)(r,  (o)  denotes  the  far-zone  spectrum  whereas  S(0)(r,  co)  denotes  the  source 
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spectrum.  We  note  that  although  most  of  our  interest  is  in  the  diagonal  component  of 
the  cross-spectral  density  function,  no  propagation  calculation  can  be  performed 
without  the  explicit  knowledge  of  the  full,  complex-valued,  cross-spectral  density. 

A  quantity  closely  related  to  the  spectrum  is  the  radiant  intensity^ /(u;  a>) 
which  is  a  measure  of  the  amount  of  energy  per  unit  frequency  radiated  at  frequency 
CO  in  a  solid  angle  d£2  about  a  direction  specified  by  unit  vector  u.  It  is  defined  by 

J(u\co)=  Urn  R^WiRu,Ru;co)  .  (1.16) 

R-*— 


It  is  clear  from  the  last  two  equations  that  the  cmly  difference  between  the  spectrum  in 
the  far  zone  co)  and  the  radiant  intensity  7(ii;  to)  is  the  scaling  factor  R^.  For 
this  reason  we  will  use  both  the  spectrum  aiKl  the  radiant  intensity  when  investigating 
spectral  effects. 

The  radiated  power  at  frequency  co  is  given  by  the  formula 


/»(©)=  J  /(«;co)df2  .  (1.17) 

(4X) 

where  (4ti)  denotes  integration  over  all  solid  angles.  The  radiated  power  is  an 
important  quantity  that  we  encounter  in  calculations  of  radiation  efficiency  and 
directivity.  We  will  sometimes  refer  to  this  quantity  as  the  total  power^"^  ai  frequency 
(0. 

The  spectrum,  the  radiant  intensity  and  the  power  are  the  main  quantities  of 
interest  in  this  thesis.  Unless  otherwise  noted,  we  will  be  concerned  with  their  values 
for  observation  points  in  the  far  zone.  According  to  Eqs.  (1.15),  (1.16)  and  (1.12) 
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we  may  express  the  far-zone  spectram,  and  the  radiant  intensity  in  terms  of  the  cross- 
spectral  density  of  the  source,  viz.. 


=  W^^\-ku,ku\(o)  ,  (1.18) 

R 

J(u;Q))  =  (2iikfcos^ew^°\-ku,lcu-,(o)  .  (1.19) 

For  general  propagation  in  source-free  region,  it  follows  from  Eqs.  (1.5) 
and  (1.6)  that  the  cross-spectral  density  function  satisfies  the  double  Helmholtz 
equation 


(Vf-t-fc2)(v2-hit2)w'(rj,r2;a))  =  0  ,  (1.20) 

where  is  the  Laplacian  operator  with  respect  to  rj.  It  is  straightforward  to  show 
that  a  general  solution  of  Eq.  (1.20)  expressing  the  cross-spectral  density  of  the  field 
in  terms  of  the  cross-spectral  density  of  the  source  is  given  by 

W(r^,r2;Q))  =  JJ W^°\si,S2;(o)G\r^,s^;a)Gir2,S2;(o)d\ d\  ,  (1.21) 

where  /i  =  1, 2,  3,  according  to  the  dimensionality  of  the  problem,  and  G  is  the  ftee 
space  Green’s  function  for  the  given  boundary  conditions.  In  particular,  for 
propagation  in  a  linear  system,  we  simply  replace  in  Eq.  (1.21)  the  Green's  function 
G(r,  s;  to)  by  the  impulse  response  of  the  system. 
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Representations  and  models  of  the  cross*spectral  density 

Another  item  of  resemblance  between  the  cross>spectral  density  and  the 
mutual  coherence  function  is  that  the  cross-spectral  density  can  be  normalized  in  a 
way  similar  to  Eq.  (1.3).  The  normalized  form  is  known  as  the  complex  degree  of 
spatial  coherence  which  is  defined  by^ 


//(rj.rj;©) 


Wir,,r^\a)) 


(1.22) 


and  which  satisfies  the  constraint 

|^(rpr2;©)|si  .  (1.23) 

Here  the  lower  limit  corresponds  to  spatial  incoherence  at  fiequcncy  ©  while  the 
upper  limit  conesponds  to  complete  spatial  coherence  at  that  frequency.  It  follows 
from  Eqs.  (1.15)  and  (1.22)  that  the  cross-spectral  density  can  be  expressed  in  the 
general  form 


H^(rj,r2;©)  =  .^S(ii^.^5(r^I^/i(rj,r2;©)  .  (1.24) 

Three  special  cases  that  are  often  encountered  in  coherence  theory  are  derived  from 
Formula  (1.24):  the  cross-spectral  densities  of  the  homogeneous^^,  Schell-nxxiel^ 
and  quasi-homogeneous^^  sources. 

A  homogeneous  source  is  represented  by  the  cross-spectral  density 
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W(r^,T^\0))  =  F{r^-r{,(0)  .  (1.25) 

This  expression  implies  that  the  source  is,  in  principle,  infinitely  large.  We  point 
out  that  this  difficulty  can  be  removed  if  we  assume  that  the  function  F  vanishes 
outside  a  finite  domain  o. 

A  Schell-model  source  is  characterized  by  a  cross-spectral  density  whose 
degree  of  correlation  depends  on  the  two  spatial  arguments  only  through  their 
difference,  i.e., 


W(r^,r^\oj)  =  ^S(r^,0))^S{r^;a>)^ir2 -r^,co)  .  (1.26) 

This  form  of  the  degree  of  correlation  is  a  natural  choice  in  nx>st  practical  situations. 
When  we  describe  secondary  sources,  it  is  reasonable  to  assume  that  the 
fluctuations  of  a  source  in  a  given  neighborhood  are  a  result  of  the  radiation  from  a 
single,  spatially  coherent,  primary  radiator.^^  When  the  separation  between  the 
points  r2  and  rj  is  larger,  the  fluctuations  reflect  contributions  from  noore  than  one 
primary  radiator  and  hence  the  fields  are  only  partially  correlated.  A  similar 
argument  can  be  made  regarding  primary  sources.  Here  one  has  to  take  into 
account  cooperative  effects  that  depend  on  the  nature  of  the  radiating  system.^^ 
Supplementing  the  choice  of  the  degree  of  correlation  in  this  model,  the  choice  of 
the  spectral  intensity  S(r,(a)  as  a  function  of  position  effectively  determines  the 
physical  size  of  the  source. 

(^uasi-homogeneous  sources  are  a  very  popular  special  case  of  the  Schell- 
model  class.  They  are  used  in  many  investigations  in  coherence  theory  for  reasons 
that  become  apparent  from  their  mathematical  and  physical  properties.  Their  cross- 
spectral  density  is  characterized  by  a  sharply  peaked  degree  of  correlation  and  a 
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broad  intensity  profile  that  does  not  change  appreciably  over  a  typical  correlation 
length  (see  Fig.  1.2).  As  a  result  of  this  relation  between  the  intensity  profile  and 
the  degree  of  correlation,  the  cross-spectral  density  may  be  approximated  by  the 
following  expression 


~r^\0))  .  (1.27) 

This  approximation  may  not  appear  advantageous  until  one  employs  it  in  calculating 
far- zone  expressions  for  the  cross-spectral  density,  the  spectrum  or  the  r^'^’ant 
intensity.  Using  the  transformation 

^ 'i)A  ’  ^2  “ 

with  unity  Jacobian,  it  is  easy  to  show  that  the  spatial  Fourier  transform  of  the 
cross-spectral  density  can  be  factexed  in  the  form 


(1.29) 


In  Eq.  (1.29)  5  and  p.  are  the  spatial  Fourier  transforms  of  the  intensity  and 
correlation  respectively,  which  for  planar  sources  are  given  by 


(1.30) 
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Figure.  1.2:  The  relation  between  the  intensity  int>file  and  the  spatial  coirelaicm  of  a 
quasi-hotnogeneous  source. 

and 

=  .  (1.31) 

(2Jt) 


I 


The  fact  that  the  spatial  Fourier  transform  of  the  cross-spectral  density  of  a  quasi- 
homogeneous  source  factorizes  into  a  product  of  two  spatial  Fourier  transforms  of 
half  the  dimensions  is  a  significant  simplification  of  the  mechanics  of  the 
propagation  calculation,  but  at  this  point  we  should  also  point  out  that  a  large  class 
of  physical  sources  can  actually  be  characterized  as  quasi -homogeneous  sources. 
In  particular,  all  sorts  of  thermal  sources  have  relatively  short  spatial  correlation 
length  and  at  the  same  time  possess  an  alnx>st  uniform  intensity  profile,  which  is 
precisely  the  description  of  a  quasi-homogeneous  source. 
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The  different  model  sources  we  have  described  represent  some  of  the 
physical  properties  of  common  light  sources.  Another  description  of  partially 
coherent  sources  can  be  made  using  coherent-mode  representation,'^^  which  is 
based  on  the  expansion  of  the  cross-spectral  density  into  a  sum  of  mutually 
uncorrelated  coherent  fields  (modes).  Before  we  proceed  with  the  formulation  of 
coherent-mode  representation,  we  point  out  some  of  the  properties  of  the  cross- 
spectral  density  function  on  which  the  representation  is  based. 

The  cross-spectral  density  is  an  hermitian  quantity,  i.e.. 


W(rj,r2;6))  =  W*(r2,rj;G))  . 


(1.32) 


It  is  also  a  non-negative  definite  quantity^,  explicitly. 


r^;(0)f\r^)fir^)d''r^  d\  ^  0  ,  (1.33) 

o 

where  n  has  the  dimensionality  of  the  source  and  /  is  any  square  integrable 
function.  Under  very  general  conditions  (in  fact,  in  all  known  practical  cases)  the 
cross-spectral  density  also  satisfies  the  relation 


jj\Wir^,r^’,Q)fd\  d\  <  oo  .  (1.34) 


When  a  continuous  function  of  two  variables  (rj  and  r-^  satisfies  conditions  (1.32), 
(1.33)  and  (1.34)  it  can  be  expressed  as  an  absolutely  and  uniformly  convergent 
series  known  as  the  Mercer  expansion 
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(1.35) 


where  the  index  m  has  the  same  dimensionality  as  that  of  and  r2,  are  the  non- 
negative  eigenvalues  and  are  the  orthogonal  eigenfunctions  of  the  Fredholm 
integral  equation 


The  eigenfunctions  can  be  made  oithonomud,  i.e.. 


(1.36) 


j ^•(r;o))^.(r,6»d''r  =  6,..  ,  (1.37) 

<T 


where  ^  is  the  Kronecker  symbol. 

It  follows  from  the  definition  of  the  degree  of  correlation  Eq.  (1.22)  that  a 
corss-spectral  density  which  consists  of  a  single  eigenfunction  is  fully  coherent,  i.e. 


Infrj.rj;© 


=  1. 


(1.38) 


The  magnitude  of  the  degree  of  conelation  of  a  source  consisting  of  several  modes 
reduces  with  increasing  number  of  modes.  The  rate  of  decrease  depends  on  the 
particular  distribution  of  the  eigenvalues.^ 
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Coherent-mode  representation  which  we  discussed  earlier  seems  to  be  a 
more  complicated  formulation  than  the  closed  form  representation  of  partrially 
coherent  fields  which  we  considered  earlier.  It  is  however  important  to  note  the 
essential  differences  between  these  representations.  Whenever  both  representations 
are  available,  it  is  usually  simpler  to  handle  most  calculations  using  the  full 
functional  form  of  the  model  source.  However,  when  the  coherent-mode 
representation  is  available  and  the  source  is  relatively  spatially  coherent,  the  number 
of  eigenfunctions  in  the  sum  is  manageable  and  it  not  only  gives  a  picture  of  the 
source  in  terms  of  a  superposition  of  coherent  fields  but  it  also  facilitates  somewhat 
simpler  calculations  because  each  one  of  the  eigenfunctions  can  be  propagated 
separately  as  a  fully  coherent  field.^2  xhe  advantage  of  propagating  each 
eigenfunction  separately  is  used  in  inverse  problems,'*^  and  in  analyzing  the  spatial 
content  of  the  illumination.'^  In  chapter  5  we  consider  a  partially  coherent  source 
and  analyze  correlation  effects  on  the  spectrum  using  coherent-mode  representation. 
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Introduction 

The  concept  of  radiation  efficiency  is  used  extensively  in  various  fields  of 
electromagnetic  theory.  ^  It  is  usually  a  measure  of  the  effectiveness  of  the  producdon 
of  radiation  by  a  given  mechanism  which  is  taken  as  the  ratio  of  the  emitted  energy  to 
the  energy  used  by  the  system.  In  many  applications  the  radiation  efficiency  also  takes 
into  account  the  spadal  or  spectral  distribution  of  the  fields. 

In  this  chapter  we  consider  the  total  flux  emitted  by  a  partially  coherent  source  at 
a  fixed  frequency  co,  and  compare  it  with  the  source  integrated  intensity  at  that 
frequency.  We  begin  our  discussion  with  the  analysis  of  the  radiation  efficiency  of 
planar  Gaussian  Schell-model  sources.  This  example  gives  a  unique  insight  into  the 
relationship  between  the  size  of  the  source  intensity  profile  and  the  spatial  correlation 
length  of  the  light  fluctuations.  We  then  examine  the  radiation  efficiency  of  three- 
dimensional  primary  sources.  The  treatment  of  certain  aspects  of  quasi-homogeneous 
sources  requires  the  extension  of  the  mathematical  definition  of  the  radiation  efficiency. 
After  considering  the  radiation  efficiency  of  a  uniform  coherent  spherical  source  we 
conclude  this  chapter  by  presenting  a  method  for  choosing  the  optimal  spatial 
correlation  which  maximizes  the  radiation  efficiency  of  a  three-dimensional,  primaiy 
sources  of  a  prescribed  intensity  profile. 
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2.2  The  radiation  efficiency  of  planar  Gaussian  Schell-model 
sources 

The  radiation  efficiency  of  partially  coherent  sources  was  studied  by  several 
investigators.^  Initial  work  focused  on  quasi-homogeneous  planar  sources.  Among 
planar  sources,  the  properties  of  Gaussian  Schell-model  sources  have  been  studied 
most  extensively .3  This  fact  may  be  attributed  to  two  main  reasons:  physically  this 
class  of  sources  is  closely  related  to  the  radiation  from  certain  types  of  lasers^  and 
mathematically  the  Gaussian  intensity  profile  and  Gaussian  correlation  functions  are 
convenient  forms  in  the  analysis. 

Wc  define  the  radiation  efficiency  of  secondary  planar  sources  by  the  formula 


C(©)  = 


<PiQ)) 

js(r;(i))d^r  ’ 

D 


(2.1) 


where  D  is  the  domain  occupied  by  the  source  and  0  is  the  total  emitted  flux  at  a 
frequency  O),  i.e.. 


0i(o)=  jj{u;a))dQ  .  (2.2) 

(2k) 

In  Eq.  (2.2)  J  is  the  radiant  intensity  and  (lit)  denotes  integration  over  the  half  space 
z  >  0  as  shown  in  Fig.  2.1. 

The  Gaussian  Schell-model  source  is  characterized  by  the  cross-spectral  density 
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Figure  2.1:  Illustrating  the  notation  used  in  this  chapter.  The  source  is  located  in 
the  plane  z  =  0  and  a  typical  observation  point  P  is  in  a  direction  specified  by  unit 
vector  n  and  at  a  distance  R  (from  the  origin)  assumed  to  be  in  the  far  zone. 


VV(rj,r2;6))  =  S^°\ci))cxp 


‘•A  . 


exp 


2(7? 


(2.3) 


where  Oj  is  the  rms  width  of  the  intensity  profile  and  is  the  rms  correlation  length 
and  5(^1((0)  is  the  source  spectral  distribution.  Although  the  cross-spectral  density  in 
£q.  (2.3)  implies,  in  principle,  an  infinite  source,  wc  note  that  whenever  the  maximum 
physical  extension  of  the  source  L  say,  satisfies  the  relation  L  »  Gj  then  the  Gaussian 
intensity  profile  is  essentially  the  blocking  function  which  defines  the  size  of  the 
effective  source. 

To  obtain  an  expression  for  the  radiation  eftlciency  of  the  Gaussian  Schell- 
rrKxiel  source  we  first  evaluate  the  denominator  of  Eq.  (2.1).  Using  £q.  (2.3)  with 
rj  =  r2  =  r,  we  have 
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D 


((D)j*  ddj  re\p^-r^/2a^j^dr 


=  27caj5^°^(0))  . 
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(2.4) 


Next  we  derive  an  expression  for  the  radiant  intensity  J(u;  to)  using  the  four¬ 
dimensional  spatial  Fourier  transform  of  the  cross-sepctral  density  [cf.  Eq.  (1.19)]. 
We  obtain 


where 


W(-kSj^,ks^-,0}) 


1 

(2)7 


sin  "  e/2) . 


(2.5) 


‘  ^  4.  * 

7' 


(2.6) 


On  substituting  Eq.  (2.5)  in  Eq.  (M9)  we  find  that  the  radiant  intensity  generated  by 
the  source  is  given  by 


/(«;  CO)  =  S^^\Q))k^<J^a^  cos^6  exp(-x^a^  sin^  d/lj  .  (2.7) 

The  expression  for  the  total  flux  is  obtained  by  substituting  Eq.  (2.7)  into  Eq.  (2.2),  it 
is  then  given  by  the  formula 


<I>iO))  =  S^^\co)[kcaif  Jcos^0  exp(-i^<7^sin^0/2)df2  .  (2.8) 

(2ii) 
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To  perform  the  integration  we  define 


cos9  =  x 
ka/yl2  =  ^  . 


(2.9) 


Using  Eqs.  (2.9)  in  Eq.  (2.8)  we  obtain 


<D((o)  =  27t5^®'^(co)(toO;f  jx^exp[-^^(l  -x^)]dx  .  (2.10) 


which  can  be  simplified  to  the  form 


0(0))  =  2nS^^\a))af 


1- 


(2.11) 


On  substituting  Eqs.  (2.4)  and  (2.11)  into  Eq.  (2.1)  we  finally  obtain  the  following 
expression  for  the  radiation  efficiency  of  planar  Gaussian  Schell-model  sources: 


C(©)  =  l-D(^)/^  ,  (2.12) 


where 


D({)  =  exp(-42)Jexp(i2)dr 

0  (2.13) 


is  the  Dawson  integral^  (see  Fig.  2.2). 
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Figure  2.2:  Graphical  rqiresentation  of  D(^/^  where  D(i^  is  the  E>awson  integral. 


The  radiation  efficiency  given  by  Eq.  (2.12)  is  a  function  of  the  source  intensity 
profile  and  its  spatial  correlation  length  through  Eqs.  (2.6)  and  (2.9).  In  particular,  we 
note  from  Eq.  (2.6)  that  there  are  classes  of  different  Gaussian  Schell-model  sources 
that  have  the  same  radiation  efficiency.  Specifically,  the  characteristic  quantity  o  which 
determines  the  radiation  efficiency  is  a  sum  of  a  term  corresponding  to  the  source  rms 
intensity  and  a  term  corresponding  to  the  correlation  length  of  the  source  fluctuations. 
It  is  clear  from  this  equation  that  there  are  many  combinations  of  different  correlation 
lengths  and  different  intensity  profiles  that  lead  to  the  same  value  of  7  and  hence  to  the 
same  radiation  efficiency.  These  equivalent  classes  of  partially  coherent  sources  were 
first  discussed  in  connection  with  the  spatial  distribution  of  the  radiation  produced  by 
Gaussian  Schell-model  sources.^  The  radiation  efficiency  in  Eq.  (2.12)  is  independent 
of  the  spectrum.  Its  only  dependence  on  frequency  is  due  to  the  wavenumber  k  and  the 
possible  frequency  dependence  of  the  parameters  and  C].  In  this  chapter  we  take 
these  parameters  to  be  independent  of  frequency. 
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Figure  2  J:  The  radiation  eHlciency  C((o)  as  a  function  of  the  normalized  spatial  correlation 
length  kxsc.  The  three  curves  indicate  that  the  radiation  efficiency  increases  rapidly  with  an  increase 
in  the  source  rms  width. 

Before  we  examine  the  equivalent  classes  of  planar  Gaussian  Schell-model 
radiators  we  consider  separately  the  effects  of  the  source  size  (rms  intensity)  and  the 
spatial  correlation  length  on  the  radiation  efficiency.  In  Fig.  2.3  we  show  curves 
representing  the  radiation  efficiency  of  three  small  sources.  The  intensity  profiles  are 
characterized  by  rms  widths  Oj  =  1/it,  2/it,  and  10//t  respectively.  We  note  that  for 
each  one  of  the  sources  the  radiation  efficiency  increases  with  increasing  correlation 
length.  The  asymptotic  value  of  the  radiation  efficiency  in  each  case  may  be  obtained 
directly  from  Eqs.  (2.6)  and  (2.9),  by  setting 

^2->2(ittT/f  .  (2.14) 

In  a  similar  manner  we  show  in  Fig.  2.4  the  radiation  efficiency  as  a  function  of 
the  source  intensity  profile  (or  the  effective  source  size).  We  note  that  the  radiation 
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Figure  2.4:  The  radiation  efficient  C(a))  as  a  function  of  the  width  of  the  source  intensity 
profile  kxfi.  The  three  curves  indicate  that  the  radiation  efficiency  increases  rapidly  with  an  increase 
in  the  source  correlation  length. 

efficiency  increases  rapidly  with  increasing  width  of  the  intensity  profile  and  effectively 
reaching  its  asymptotic  value  for  a  relatively  small  source. 

The  results  shown  in  the  last  two  figures  are  combined  in  Fig.  2.5  where  the 
radiation  efficiency  is  shown  as  a  function  of  both  parameters  koi  and  kOc.  In  this 
picture  an  equivalent  class  of  Gaussian  Schell-model  sources  is  represented  by  the 
infinite  number  of  sources  specffied  by  pairs  of  values  of  kOi  and  kc^  that  lie  on  the 
intersection  of  any  plane  parallel  to  the  kaj-ka^.  plane  and  the  plotted  radiation 
efficiency  surface.  Several  such  curves  are  shown  in  Fig.  2.6. 

Before  we  conclude  the  discussion  of  the  radiation  efficiency  of  planar  sources 
we  point  out  a  few  physical  considerations.  Most  secondary  sources  of  practical 
interest  are  likely  to  have  spatial  dimensicms  that  are  much  larger  than  a  wavelength.^  If 
we  denote  a  typical  source  dimension  by  L,  we  then  have  L»  X,  and  according  to 
Fig.  2.4  the  radiation  efficiency  is  nearly  independent  of  the  intensity  profile  for  this 
range  of  rms  source  intensity.  In  addition,  if  we  note  that  the  smallest  natural  spatial 
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Figure  2.5:  The  radiation  efficiency  of  Gaussian  Schell-model  sources  as  a  function  of 
the  rms  intensity  and  the  conslation  length. 


correlation  length*  is  on  the  order  of  Xfl  corresponding  to  ka^  >  3,  then  it  follows 
ti-om  Fig.  2.3  that  the  radiation  efficiency  of  most  planar  Gaussian  Schell-model 
sources  is  over  80%. 


Figure  2.6:  Contours  of  constant  radiation  efficiency  of  Gaussian  Schell-model  sources. 
Each  curve  corresponds  to  a  class  of  different  partially  coherent  sources  that  have  the  same 
radiation  efficiency. 
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2.3  Optimal  radiation  efficiency  of  three-dimensional 
primary  sources 

In  this  section  we  analyze  the  radiation  produced  by  arbitrary  three-dimensional 
scalar  source  distributions,  whose  degree  of  correlation  at  any  pair  of  points  depends 
only  on  the  (oriented)  separation  between  the  points.  Working  in  a  Fourier  transform 
space,  we  establish  an  operational  method  for  calculating  the  spatial  correlations  that 
maximize  the  radiation  efficiency  of  such  sources.  The  optimal  correlation  function 
which  we  derive  is,  in  general,  not  unique,  and  it  is  found  to  depend  on  the  distribution 
of  the  optical  intensity  across  the  source. 

We  begin  the  analysis  by  generalizing  the  definition  of  the  radiation  efficiency 
Eq.  (2.1)  to  the  form 


C(a))  = 


0(a>) 


where 


N  =  JSg(r)dV. 

D 


Here  D  is  the  source  domain  and  17  is  given  by 


(2.15) 


(2.16) 


27  = 


(2.17) 
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The  use  of  the  effective  volume  iZ),  rather  than  the  actual  volume  D,  enables  us  to 
discuss  infinite  sources  whose  optical  intensity  distribution  falls  off  sufficiently 
rapidly.^  The  advantage  of  this  formulation  will  become  apparent  later.  In  Eqs.  (2.16) 
and  (2.17)  the  subscript  Q  represents  a  primary  source  distribution. 

Another  difference  between  the  calculations  of  the  radiation  efficiency  of  planar 
and  three-dimensional  sources  is  in  the  evaluation  of  the  total  flux  which,  in  the  present 
case,  is  given  by  the  formula 


<r>(a))=  Jy(«;o))df2  ,  (2.18) 

(4k) 

where  (4«)  denotes  integration  over  all  solid  angles. 

The  radiant  intensity  generated  by  three-dimensional  primary  sources  of  any 
state  of  coherence  is  given  by  Eq.  (3.9)  of  Ref  1 1  (a). 

/(«;©)  =  (2nfW{-ku,ku;G>),  (2.19) 

where  W  is  the  six-dimensional  spatial  Fourier  transform  of  the  source  cross-spectral 
density  i.e.. 


Jj W'e(r,,r2;6))exp[-i(/j  •  Tj  -h  /j  •  r2)]d\  d\ .  (2.20) 
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The  definition  of  the  radiation  efficiency  in  Eq.  (2.15)  is  consistent  with  that 
given  in  Eq.  (2. 1)  in  that  the  efficiency  of  primary  sources  is  also  bounded  by 

O^C(0))^1.  (2.21) 

This  relation  can  be  derived  by  considering  any  source  cross-spectral  density  of  the 
form 


WQ(r^,r2;o))  = 


(2.22) 


Since  the  Hearee  correlation  is  normolLzed  so  that 


(2.23) 


it  follows  from  Eqs.(2.17),  (2.19)  and  (2.20)  that 


ly(«;to)|  <  JJ^Sg(rj;©)^Sg(r2;fi))  d\  d\  =  'DN  .  (2.24) 

The  normalization  of  the  radiation  efficiency  Eq.  (2.21)  can  now  be  obtained  on 
substituting  Eq.  (2.24)  in  Eqs.  (2.18)  and  (2.15). 

The  result  expressed  by  Eq.  (2.21)  is  valid  for  all  primary  sources  regardless  of 
their  state  of  coherence  or  the  distribution  of  their  optical  intensity.  It  is  clear  that  the 
lower  limit  C(tD)  0  is  obtained,  for  example,  in  the  limit  of  complete  spatial 
incoherence*®  since,  according  to  Eqs.  (2.19)-(2.20),  for  such  a  source  the  radiant 
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Figure  2.7:  Normalized  radiant  intensity  produced  by  a  uniform,  co-phasal,  fiilly  coherent 
^jheiical  source  of  radius  a. 

intensity  /(ii;  to)  vanishes.  Other  examples  of  this  type  are  certain  fully  coherent  and 
partially  coherent  non-radiating  source  distributions.^^ 

We  note  that  lower  limit  of  the  radiation  efficiency  for  three-dimensional 
primary  sources  occurs,  as  in  the  case  of  secondary  planar  sources,  in  the  limit  of  very 
short  spatial  correlation  lengths.  However,  as  we  show  in  ihc  following  example, 
unlike  secondary  planar  sources,  the  radiation  efficiency  is  not  always  maximized  in  the 
limit  of  long  correlation  length. 

Consider  a  uniform,  co-phasal,  fully  coherent  spherical  source  of  radius  a.  The 
radiant  intensity  produced  by  such  a  source  is  given  by  [Ref.  1 1(a)  Eq.  (4.15)] 


J(u-,(a)  =  (4rca^/3f  Sg(©)[3yi(to)//taf  ,  (2.25) 

see  Fig.  2.7.  In  Eq.  (2.25)  Jiix)  is  the  spherical  Bessel  function  of  the  first  order  i.e.. 
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Jl(x)  = 


SillX 


cosx 

X 


(2.26) 


In  the  limit  of  small  argument  x  «  1  we  have  Jiix)  *  x/3.  It  then  follows  from  Eq. 
(2.25)  that  the  radiant  intensity  generated  by  a  spherical  source  of  this  type,  whose 
radius  a  is  smaller  than  the  wavelength  is 


7(ii;m)  =  (4na^/3fsQ(r,0))  =  D\{r\(o)  ,  .  (2.27) 


where  D  is  the  volume  of  the  source.  On  the  other  hand,  for  a  source  of  this  type 
N  =  DSq(©)  and  D  =  D.  Then  according  to  Eq.  (2.15),  the  radiation  efficiency 
C((D)  1  as  fl  -4  0.  This  result  shows  that  the  maximum  radiation  efficiency  is 

achieved  by  a  uniform  coherent  source  in  the  limit  as  its  linear  dimensions  are  small 
compared  with  the  wavelength.  When  the  dimensions  of  the  source  increase  the 
radiation  efficiency  decreases  and  actually  vanishes  for  radii  a  such  that  ila  is  a  zero  of 
the  Bessel  function  This  result  is  in  contrast  with  the  result  for  the  radiation 
efficiency  of  planar  secondary  sources  that  achieve  maximum  efficiency  for  fully 
coherent  field  distributions  and  large  spatial  dimensions  (see  Fig.  2.3). 

We  now  wish  to  find  the  spatial  correlation  function  gQ(rj  -  r^,  ©)  which 
maximizes  the  radiation  efficiency  C((i})  for  a  given  primary  source  with  an  arbitrary 
(but  known)  intensity  distribution  SQ(r;  co)  and  domain  D.  Since  the  quantities  N  and 
D  are  then  fixed,  we  may  confine  ourselves,  in  view  of  Eq.  (2.15),  merely  to  the 
optimization  of  gQ(r]  -  r^,  co)  in  such  a  way  that  the  total  flux  <I>((0)  radiated  by  the 
source  attains  its  absolute  maximum  value. 

According  to  F.qs,  (2.18)-(2.20)  the  radiated  flux  at  frequency  ©  is  given  by  the 
expression 
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In  Eqs.  (2.30)  and  (2.31)  the  spatial  integrations  are  performed  over  the  domains  £>i 
and  D2  which  are  determined  from  the  original  source  domain  D  through  the  coordinate 
transformation  Eq.  (2.29).  For  example,  it  is  evident  from  Eq.  (2.28)  that  the  modulus 
of  r'  does  not  exceed  the  maximum  linear  dimension  of  D.  The  function  /?(r’;  (o)  is 
related  to  the  source-averaged  correlation  function  employed  in  earlier  studies  of 
radiation  by  partially  coherent  planar  sources.*'* 
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The  next  step  in  the  analysis  is  the  integration  of  Eq.  (2.31)  with  respect  to  Cl. 
Since  the  integral 


/  = 


J  exp(-ifeu  •  r')dCl 
(4k) 


(2.32) 


has  no  preferred  direction  in  space,  we  can  rotate  the  coordinate  system  so  that  r'  is  in 
the  direction  of  the  polar  axis  of  a  spherical  polar  coordinate  system,  and  the  integral  is 
readily  performed. 


J  cxp{-iku  ■  r')dC2 


2k  k 

=  J d^jexp(-i/l:r'cos0)  sinddd 
0  0 


-1 


(2.33) 


=  4k 


sin  to-’ 
kr' 


where  r'  =  Ir"!.  On  substituting  Eq.  (2.33)  into  Eq.  (2.31)  we  find  that 


<*>(<y)=  jG(r’)gg(r’)dV  , 
^2 


(2.34) 


in  which  we  have  defined  the  auxiliary  function  Qir*;  to)  by 


Q(r';o))  =  47t/?(r’;0)) 


sinkr' 
kr'  ' 


(2.35) 
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Here  the  functions  /?(r';  to)  and  Q(r';  (O)  depend  explicitly  on  the  frequency  (o,  but  we 
will  suppress  this  dependence  in  the  rest  of  this  chapter. 

Since  the  source-averaged  intensity  function  /?(r')  vanishes  outside  the  domain 
D2,  the  same  holds  true  for  Q(r')  and  the  integration  in  Eq.  (2.34)  can  be  formally 
taken  over  the  infinite  three-dimensional  r'-space. 

We  now  introduce  the  three-dimensional  spatial  Fourier  transforms  of  Q(r’)  and 
^q(^)  by  the  formulas 


Q(K)  =  j Q(r')exp(-iK ■  r')d^r'  .  (2.36) 

gQ(f^)  =  -^^jgQ(r')cxp{-iK  r')d^r'  .  (2.37) 

On  invoking  Parseval’s  theorem'^  in  Eq.  (2.34),  we  noay  then  express  the  total  radiated 
flux  by** 


0(©)  =  (2n)^ J Q{K)gQ{K)d^K  .  (2.38) 

At  this  stage  we  focus  our  attention  on  the  mathematical  properties  of  the  functions  Q 
and  gQ.  Because  the  function  Qif)  is  real  and  inversion  symmetric  with  respect  to  the 
origin,  we  can  establish  the  relations 

Q\K)  =  Q{-K)  =  Q{K)  .  (2.39) 
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It  follows  from  Eq.  (2.39)  that  Q  is  also  real  and  inversion  symmetric  function  of  K 
with  respect  to  the  origin  in  X-space.  The  other  factor  in  the  integrand  of  Eq.  (2.38), 
gg,  is  the  Fourier  transform  of  the  complex  degree  of  spatial  coherence  of  the  source 
fluctuations.  According  to  Bochner’s  theorem^^,  a  necessary  and  sufficient  condition 
for  gg(r')  to  be  a  correlation  coefficient^*  is  that  it  is  non-negative  definite  and  it  is 
equal  to  unity  at  the  origin.  In  the  Fourier  domain  these  requirements  are  equivalent  to 


0, 


(2.40) 


and 


gQiO)  =  \gQ{K)d^K  =  \  .  (2.41) 

It  follows  from  Eq.  (2.38)  and  from  the  properties  of  the  two  factors  in  the  integrand  of 
that  equation,  that  there  is  a  simple  operational  method  for  finding  a  spatial  correlation 
function  go(r')  which  maximizes  the  total  radiated  flux  d>{co)  and  consequently  also 
maximizes  the  radiation  efficiency  C(a)).  In  its  simplest  form,  the  basic  principle  of 
this  method  is  to  determine  the  point  A  =  Hlq  at  which  Q{K)  assumes  its  absolute 
maximum  value;  the  optimal  correlation  function  then  corresponds  to 

g{K)  =  5^^\K-KQ)  .  (2.42) 

where  5^^^  is  the  three-dimensional  Dirac  delta  function.  This  is  only  a  particular 
solution  of  the  maximization  problem  which  has  several  interesting  consequences. 
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The  first  consequence  of  this  solution  is  that  an  upper  bound  on  the  amount  of 
flux  that  a  partially  coherent  source  can  generate  is  given  by 

=  (2.43) 

The  radiation  efficiency  corresoonding  to  this  maximum  flux  is  readily  obtained  by 
substituting  Eq.(2.43)  into  Eq.  (2.15). 

The  second  consequence  of  Eq.  (2.42)  is  that  at  least  one  solution  to  the 
maximization  problem  is  a  fully  coherent  source.  This  follows  from  taking  the  inverse 
Fourier  transform  of  Eq.  (2.42),  i.e., 

gg(r‘)=:exp(iA:or').  (2.44) 

It  is  important  to  note  at  this  stage  that  the  optimal  correlation  indeed  corresponds  to  a 
coherent  source,  but  the  phase  of  the  correlation  is  not  uniform  as  can  be  seen  from  Eq. 
(2.44).  The  linear  phase  factor  in  this  equation  depends  on  the  optical  intensity 
distribudon  5g(r)  through  the  maximum  point  Kq  of  the  function  Q- 

The  third  consequence  of  our  solution  to  the  maximization  problem  is  that  the 
optimal  spatial  correlation  is,  in  general,  not  unique.  Different  source  correlation 
coefficients  gg(r’)  may  yield  the  same  maximum  radiation  efficiency.  Moreover,  as  we 
pointed  out  earlier,  Q  is  inversion  symmetric  and  hence  the  point  K  =  -Kq 
corresponds  to  a  maximum  whenever  K  =  Kq  does.  Consequently  5^^\K  +  Kq)  is 
also  a  particular  solution  of  the  maximization  problem.  Either  one  of  the  two 
unimodular  source  correlation  functions  gQ(r')  can  then  be  regarded  as  the  optimal 
complex  degree  of  spatial  coherence.  Similarly,  a  linear  combination  of  the  form 
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g{K)  =  -  Ko)  +  q8^'^\K  +  K^)  ,  (2.45) 

where  p  +  q  =  1,  is  also  in  agreement  with  Eqs.  (2.40)  and  (2.41)  and  produces, 
according  to  Eq.  (2.38),  the  maximum  possible  radiant  flux  <I>(o))mtx-  This  choice  of 
giK)  however  does  not  correspond  to  spatially  fully  coherent  source  distribution.  In 
fact,  Eq.  (2.45)  corresponds  to  a  degree  of  correlation  of  the  form 

gQ (r')  =  p exp[iKQ  •  r’)  +  ^exp(-iA:o  -  r')  ,  (2.46) 

whose  magnitude  is  given  by^^ 

jpexp(tArQ  •  r’)  +  <7exp(-iA'o  •  r’)j  =  +  2pqcosi2kr')  .  (2.47) 

In  general,  Q(K)  may  take  on  its  maximum  value  over  an  extended  domain 
such  as  a  surface  in  the  three-dimensional  K-space.  In  such  a  case  the  maximizing 
g(K)  may  be  distributed  arbitrarily  throughout  the  domain  provided  only  that  the 
conditions  of  Eqs.  (2.40)  and  (2.41)  are  met  Hence  an  infinite  number  of  optimal 
correlations  with  varying  degrees  of  spatial  coherence  can  then  be  specified.  In  the 
following  paragraphs  we  consider  several  physically  in^iortant  applications  where  these 
situations  occur. 

The  first  case  that  we  consider  is  that  of  quasi-homogeneous  sources.  The 
cross- spectral  density  of  such  sources  may  be  expressed  in  the  form 

W^(ri,r2;a))  =  S[(r,  +r2)/2;®)]gg(r,  -rj)  .  (2.48) 
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Wc  note  that  for  quasi-homogeneous  sources  we  may  effectively  set  Dj  *  D  and 

VS(r  +  ry2; (0)^Sir  -  ry2; Q))  =  S(r,  to)  ,  (2.49) 


in  Eq.  (2.30).  The  source-averaged  quantity  /?(r')  then  assumes  a  constant  value 


/?(r')  =  j5g(r)dV  (2.50) 

D 

It  follows  from  Eqs.  (2.50)  and  (2.35),  that  the  function  Q(r')  is,  in  this  case, 
independent  of  the  source  domain  D;  and  it  is  given  by  the  spherically  symmetric 
(isotropic)  formula 


Q{r’)  =  4nN 


sin  kr' 
kr' 


(2.51) 


We  note,  however,  that  strictly  speaking  there  should  be  a  finite  (generally  non¬ 
isotropic)  cut-off  in  the  values  of  r"  in  Eqs.  (2.50)  and  (2.51)  due  to  the  domain  D2.  In 
the  quasi-homogeneous  approximation  all  contributions  to  the  radiant  flux  <P(G))  arise 
from  a  small  neighborhood  c  O2  around  the  origin  r'  =  0  in  which  gQ(r’)  differs 
appreciably  from  zero.  We  therefore  let  the  function  Q(r')  extend  over  the  entire  space 
without  altering  the  situation  physically. 

Allowing  for  an  infinite  domain  (D2  —*  “>)  the  spatial  Fourier  transform  of 
!2(/')  is  found  by  substituting  from  Eq.  (2.51)  into  Eq.  (2.36).  One  then  finds  that 


Q(K)  =  p[5(|^:|  -k)-  (5(|/(:|  +  k)].  (2.52) 
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The  second  terai  in  Eq.  (2.52)  makes  no  contribution  to  the  flux  O(o)).  We  therefore 
set 


Q(K)  =  ^S{K-k)  ,  (2.53) 

where  =  llfl.  We  observe  that  for  quasi'homogeneous  sources  Q(K)  is  isotropic 
and  vanishes  identically  except  on  the  energy  shell  AT  =  in  the  three-dimensional  K- 
space.  Hence,  in  view  of  Eq.  (2.38),  the  quasi-homogeneous  source  will  be  totally 
non-radiating  if  its  complex  degree  of  spatial  coherence  is  such  that  g(K)  =  0  when 
K  =  k.  This  result  is  in  agreement  with  earlier  investigations  [  see,  for  example,  Eq. 
(3.11)  of  Ref.  11(b)]. 

When  the  finite  size  of  the  quasi-homogeneous  source  is  taken  into  account,  it  is 
evident  that  Q(K)  will  no  longer  be  simply  proportional  to  a  5-function  as  in  Eq. 
(2.53).  Instead  it  will  be  a  sharply  peaked,  generally  non-isotropic  function,  centered 
around  the  shell  K  =  k.  In  such  a  case  the  optimal  correlations  can  be  found  by  using 
the  general  method  we  described  above.  Nevertheless,  we  show  that  physically 
interesting  results  are  readily  obtained  with  the  help  of  Eq.  (2.53)  for  limitingly  large 
quasi-homogeneous  sources. 

Since  the  function  Q(i(()  in  Eq.  (2.53)  is  spherically  symmetric,  the  values  of 
the  Fourier  transform  ggi/^)  on  the  shell  K  =  k  can  be  chosen  in  a  number  of  ways. 
Specifically,  in  view  of  Eqs.  (2.41)  and  (2.53),  one  may  set  as  the  optimal  correlation 
fiinction 


gQiK)  =  ^F{u)5{K-k)  , 


(2.54) 
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where  K  =  Ku  and  F(m)  is  an  arbitrary  non-negative  function  that  satisfies  the 
condition 


jF(u)dn  =  l.  (2.55) 

(4ii) 

Using  the  inverse  of  Eq.  (237),  it  follows  from  Eq.  (2.54)  that  the  spatial  correlation  is 
given  by 


g(r')=  jF(u)cxp{ikur')d£2.  (2.56) 

(4x) 

In  the  case  of  an  isotropic  correlation  F(u)  =  l/47i  and  Eq.  (2.56)  yields  [cf.  Eq. 
(2.33)]  the  spatial  correlation 


gir’)  = 


sin  kr' 
kr' 


(2.57) 


This  is  an  interesting  result.  It  implies  that  in  the  limit  as  the  size  of  a  quasi- 
homogeneous  source  approaches  infinity,  the  optimal  correlations  that  are  isotropic  and 
that  maximize  the  radiation  efficiency  of  the  source  are  given  by  the  universal  function 
sin(j:)/j:  with  x  =  kr'.  Since  the  correlation  distance  associated  with  the  degree  of 
coherence  given  by  Eq.  (2.57)  is  about  X/2,  the  efficient  radiation  characteristics  of 
such  a  source  can  be  physically  understood  on  the  basis  of  an  interference  model.^o 
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Figure  2.8:  Physical  model  explaining  the  reduced  radiation  efflciency  of 

coirelated  piimary  sources  as  a  result  of  destnictive  interference  [cf.  Ref.  S(d)]. 

Consider  a  three  dimensional  source  distribution  in  which  several  small 
radiating  volume  elements  are  equally  spaced  on  an  arbitrary  axis  (see  Fig.  2.8).  When 
all  the  volume  elements  are  radiating  coherently  and  in  phase,  it  follows  that  the 
radiation  from  any  pair  of  elements  that  are  separated  by  a  distance  r»f  XP.  cancels 
everywhere  on  the  axis.  If,  however,  the  spatial  correlation  of  the  source  fluctuations 
has  an  effective  length  shorter  than  A/2,  the  destructive  interference  effects  do  not  take 
place. 

The  result  given  by  Eq.  (2.57)  is  also  closely  related  to  the  radiation  produced 
by  a  source  that  is  in  thermal  equilibrium  with  its  surroundings,  i.e.,  a  blackbody 
radiator.  It  was  shown  recently ,2*  that  a  fluctuating  source  with  this  form  of  spatial 
correlations  gives  rise  within  the  source  medium  to  a  field  whose  complex  degree  of 
spatial  coherence  is  also  given  by  the  function  sinikr')lkr'.  This,  in  turn,  is  the 
correlation  function  that  is  known  to  be  associated  with  blackbody  radiation  fields.22 
Hence  our  results  are  consistent  with  the  notion  that  a  large  uniform  source  in  thermal 
equilibrium  has  isotropic  spatial  correlations  that  lead  to  maximum  radiation  efficiency. 

As  another  illustration  of  the  operational  method  for  finding  the  spatial 
correlation  function  which  maximizes  the  radiation  efficiency,  we  consider  three 
spherically  symmetric  Schell-model  sources,  aU  of  which  have  identical  Gaussian 
optical  intensity  distributions: 
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Sg(r;®)  =  /(£0)exp(-rY2Og)  .  (2.58) 

Here  /(co)  characterizes  the  spectral  content  of  the  source  and  <7^  is  a  positive  constant. 

We  now  compare  three  types  partially  coherent  sources  having  this  intensity 
profile.  The  sources  are:  (a)  fully  coherent,  (b)  Gaussian  correlated  (optimal  correlation 
length  as  in  Ref.  2(d),  and  (c)  optimally  correlated  as  discussed  in  this  section. 

For  the  fully  coherent  source  we  take  gQir')  =  1,  whereas  the  Gaussian 
correlated  source  is  specified  by 

gg(r')  =  exp(-r'Y2o^)  ,  (2.59) 


and  the  optimal  correlation  length  is  given  by 


(2.60) 


when 


kao<—- 
2  2 


(2.61) 


The  general  optimal  correlations  can  be  found  by  determining  the  points  where 
the  spatial  Fourier  transform  of  the  function  Q(r')  corresponding  to  Eq.  (2.58)  attains 
its  absolute  maximum  value.  The  first  step  in  the  procedure  is  to  obtain  an  expression 
for  the  auxiliary  function  7?(r')  by  substituting  Eq.  (2.58)  into  Eq.  (2.30).  After 
performing  the  indicated  integration  we  obtain  the  formula 


Radiation  efficiency  and  directivity 


Chapter  2 


58 


(2.61) 


Next  we  substitute  Eq.  (2.61)  in  Eqs.  (2.35)  and  (2.36)  to  obtain  the  spatial  Fourier 
transform  QiK) : 


Q(K)  =  47c/(o»^{exp[-2(JI:  -  Kf  a^]  -  exp[-2(it  +  Kf  o^]|  .  (2.62) 


For  kCQ  <  ■\/3/2  this  function  assumes  its  maximum  value  at  the  origin  K  =  0.  When 
k<TQ  ^  V3/2,  the  radius  K  corresponding  to  the  maximum  of  Q{K)  for  any  given  Gq 
can  be  readily  found  using  numerical  methods.  We  will  denote  this  radius  by 

Making  use  of  Eqs.  (2.17)  and  (2.58)  we  fmd  that  for  all  the  three  model 
sources  that  we  arc  considering 

2W  =  iAufl{a)a^Q  .  (2.63) 

The  total  radiant  flux  generated  by  a  Gaussian  Schell-model  source  is  given  by^ 

d>(a))  =  47i(27iaga)^/((o)exp[-(k(T)^/2]  ,  (2.64) 


where 


1  _  1  1 


(2.65) 
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and  Og  is  to  be  chosen  as  described  above  [Eqs.  (2.60)  and  (2.61)]  for  maximum 
radiated  power.  In  the  case  of  the  optimally  correlated  source,  the  radiated  flux  is  given 
by  Eq.  (2.43), 

0((o)  =  i2nfQiKf^)  .  (2.66) 

Because  of  the  isotropy  of  i2(i^)  this  case,  the  correlations  that  lead  to  this  radiant 
flux  can,  in  general,  be  chosen  in  many  different  ways.  However,  when  jfctjg  <  V3/2, 
the  value  Kq=Q  is  the  unique  solution.  This  leads  to  gQ(K)  =  6  (K)  and 
consequently  the  source  is  fully  coherent  and  co-phasal,  i.e.,  gg(r')  =  l.  On 
combining  Eqs.  (2.62)  and  (2.66)  and  taking  the  limit  as  0,  we  indeed  see  that 
the  resulting  radiant  flux  <P(o))  is  identical  to  that  given  by  Eq.  (2.64)  with  k(Tg  «». 
Here  we  have  used  the  fact  that 


-exp  [-2o^(ifc  +  /02]}  =  4it(Te  .  (2.67) 


Making  use  of  the  above  results  we  may  now  calculate  the  radiation  efficiency 
defined  by  Eq.  (2.15).  For  the  fully  coherent  and  optimally  Gaussian  correlated 
sources  we  find  from  Eqs.  (2.63)  and  (2.64)  the  formula 
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Figure  2.9:  The  radiation  efficiency  of  the  three  sources.  Note  that  the  radiation  efficiency  of  the 
sources  decreases  significantly  with  increasing  source  size,  and  that  the  three  sources  have  identical 
efficiency  for  values  of  koi  ^  0.86. 


where  a  is  given  by  Eq.  (2.64).  The  radiation  efficiency  of  the  coherent  source  is 
obtained  in  the  linait  ka^  ->  oo,  and  it  is  given  by 


C(a})  =  exp 


(2.69) 


Similarly,  for  the  radiation  efficiency  of  an  optimally  correlated  source  we  obtain  firom 
Eqs.  (2.62),  (2.63)  and  (2.65)  the  expression 


C((0)  =  -  *^of  ]  -  “p[-2oe(*  +  ]}  .  (2 

where  Kq  is  the  location  of  the  maximum  of  QiK)  which  is  determined  by  numerical 
evaluation.  These  results  are  plotted  in  Figs.  2.9  and  2.10  as  a  funed^n  of  ikOg.  As 
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Figure  2.10;  Logarithmic  plot  of  the  radiabon  efficiency  of  the  three  sources. 

shown  in  Rg.  2.9,  the  three  curves  are  identical  in  the  region  ifcOg  <  V3/2,  and  take  on 
the  maximum  value  of  unity  as  kOq  — >  0.  Figure  2.10  illustrates  the  differences 
between  the  radiation  efficiencies  of  the  three  model  sources  in  the  region  kOg  >  V3/2. 

Taking  the  correlations  that  maximize  the  radiation  efficiency  to  be  isotropic, 
i.e.,  choosing 


gQiK)  =  SiK  -  ,  (2.70) 

we  readily  find  on  taking  the  inverse  Fourier  transform  of  Eq.  (2.70)  that  the  optimal 
degree  of  spatial  coherence  is  given  by 


%(r’)  = 


sin  K^r' 

v 


(2.71) 
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As  we  pointed  out,  for  small  sources  that  satisfy  0  <  kOQ  <  we  have 
f^Q  =  0  and  Eq.  (2.71)  corresponds  to  full  coherence.  On  the  other  hand,  for 
sufficiently  large  sources,  we  see  from  Eq.  (2.62)  that  the  maximum  radiation 
efficiency  is  attained  for  Kq">  k  (but  always  Kq  <  k).  Hence  in  the  limit  as 
taa->  oo,  the  optimal  correlations  given  by  Eq.  (2.71)  reduce  to  the  universal  result 
[Eq.  (2.57)],  derived  in  the  context  of  quasi-homogeneous  sources. 

2.4  Summary 

In  this  chapter  we  discussed  the  radiation  efficiency  of  secondary  and  primary 
partially  coherent  sources.  We  examined  the  dependence  of  the  radiation  efficiency  and 
the  radiated  flux  on  the  spatial  correlations  of  the  source.  We  showed  that  typical 
planar  secondary  Gaussian  Schell  sources  have  high  radiation  efficiency  which 
increases  with  an  increase  in  the  spatial  correlation  length  of  the  source  fluctuations. 
On  the  other  hand  we  demonstrated  that  three-dimensional  primary  coherent  sources  are 
not  necessarily  characterized  by  high  radiation  efficiency.  In  particular,  an  isotropic 
source  that  is  characterized  by  a  Gaussian  intensity  profile  has  maximum  radiation 
efficiency  for  a  spatial  correlation  given  by  Eq.  (2.71). 

The  method  for  choosing  the  spatial  correiation  which  maximizes  the  radiation 
efficiency  of  three-dimensional  primary  sources  is  significant  in  two  aspects.  First  it 
shows  from  the  point  of  view  of  coherence  theory  that  the  spatial  correlation  of 
blackbody  radiators  maximize  the  radiation  efficiency  of  the  source.  The  secoiKl  aspect 
has  to  do  with  the  possibility  of  producing  more  efficient  light  sources.  Most  of  the 
light  sources  currently  in  use  for  illumination  purposes  are  partially  coherent  and  hence 
even  a  small  improvement  in  their  efficiency  has  great  significance  for  the  energy 
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market  At  present  most  of  the  practical  methods  for  controlling  source  correlations 
that  we  are  aware  of  pertain  to  secondary  planar  sources,  but  new  approaches  to 
controlling  the  spatial  correlations  of  primary  sources  are  currently  being  introduced. 

Practical  methods  for  controlling  source  correlations  are  also  important  for  the 
discussion  in  the  next  chapter  were  we  consider  the  effects  of  the  source  correlations  on 
the  spectrum  of  the  radiation. 
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We  emphasize  that  the  infinite  size  of  the  source  is  only  used  for  the 
mathematical  convinience  of  extending  certain  integration  limits.  See  also  the 
discussion  following  Eq.  (2.3). 
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3.1  Introduction 

Techniques  for  modifying  of  the  spectrum  of  radiation  have  been  known  for  some  time 
both  in  optics  and  in  antenna  theory.  The  main  method  for  changing  the  spectrum  of 
field  is  to  vary  the  spectrum  of  the  source  itself,  assuming  that  the  spectrum  of  the 
radiation  is  simply  proportional  to  the  source  spectrum.  As  we  now  know,  this 
assumption  is,  in  general,  not  valid  in  the  case  of  partially  coherent  sources,  where  the 
spectrum  may  change  on  propagation. 

In  this  chapter  we  show  that  source-correlations  can  give  rise  to  interesting 
modifications  of  spectra.  We  consider  a  simple  physical  configuration,  consisting  of 
two  small  sources  which  generate  fields  of  identical  spectra,  and  we  analyze  the  effects 
of  correlation  between  the  two  sources  on  the  spectrum  of  the  emitted  radiation.  We 
show  that  spectral  lines  can  be  frequency-shifted,  made  narrower  or  broader  and  that 
several  lines  may  be  generated  from  a  single  line  by  this  mechanism.  These  results 
suggest  a  new  technique  for  modifying  spectra  in  a  desired  manner  by  controlling 
source  correlations.  In  fact,  since  these  results  were  first  published,  there  have  been 
two  reports  of  r^ectral  modulation  by  control  of  source  correlations. ^ 

The  analysis  presented  in  this  chapter  refers  to  a  simple  configuration  and  to 
observation  points  on  the  axis  of  symmetry.  The  more  general  case  of  arbitrary 
observation  point  is  discussed  in  Chapter  4. 


3.2  Radiation  from  two  small  scalar  sources^ 

Consider  light  generated  by  two  small  fluctuating  sources  located  at  points  and  P2. 
We  assume  that  the  fluctuations  are  statistically  stationary.  Let  {j2i(tD)}  and  {(22(®)} 
be  the  ensembles  that  represent  the  source  fluctuations  at  frequency  co.  Further  let 
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{ V(/*;  0))}  be  the  ensemble  that  represents  the  field  at  P,  generated  by  the  two  sources 
(see  Fig.3.1).  We  assume  that  the  spectra  of  the  two  source-distributions  are  identical 
and  we  denote  them  by  Sg(G>).  Specifically 

5e«u)  =  ((2*(6))  Q^{(o))  =  {q’^{(o)  (22(fi»)  .  (3.1) 

where  the  angular  brackets  denote  the  ensemble  average. 

The  field  produced  by  the  two  sources  at  a  point  P  is  given  by 

ikR^  ikli^ 

V(P;6))  =  (2i(6))^-K22(fi»^  .  (3.2) 

where  /?i  and  /?2  denote  the  distance  from  the  two  sources  to  the  observation  point  P 
The  spectrum  of  me  field  at  the  point  P  is  then  given  by 

Sv,(P;m)  =  Sg(<y){Pi-2  +  /2j2+29?[;ig(fi,)e^^^2-^>y/?,^]}  ,  (3.3) 

where  SR  denotes  the  real  part  and  Hq((0),  known  as  the  degree  of  spatial  coherence, 
characterizes  the  correlation  between  the  two  fluctuating  sources.  Explicitly, 

=  (G*(m)02(©))/5g(m).  (3.4) 


Spectral  modulation  with  scalar  sources 


Chapters 


Figure  3.1:  Illustrating  the  configuration.  The  two  small  sources  Qi  and  Q2  are 
place  symmetrically  about  the  axis.  A  typical  observation  point  P  is  at  a  distance  Ri 
and  R2  respectively  from  the  two  sources. 


As  we  have  explained  above  [cf.  Eq.  (1.23)],  the  degree  of  spatial  coherence  satisfies 
the  inequality 


(3.5) 

for  all  frequencies.  This  condition  is  the  main  constraint  on  the  possible  types  of 
spectral  effects  that  can  be  achieved  with  tlie  present  system. 

For  simplicity  we  will  consider  the  spectrum  of  the  emitted  radiation  at  points 
located  on  the  perpendicular  bisector,  which  we  will  refer  to  as  the  axis,  of  the  line 
joining  the  two  sources.  In  this  case  R2=  Ri  (  =  R  say)  and  Eq.  (3.3)  reduces  to^ 

Sv,(P;to)  =  -^Sg(tu)[l  +  9I;Zg(a))]  .  (3.6) 


We  note  that  when  5K/ig(o))  is  independent  of  (O,  the  spectrum  of  the  field  at  all  axial 
points  will  be  proportional  to  the  source  spectrum  Sq^ud).  This  includes  the  case  where 
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the  two  sources  are  mutually  incoherent  at  each  fiequency  IPq((o)  =  0].  It  also  includes 
the  case  when  flf/tgCa))  =  1,  which  corresponds  to  two  sources  that  are  mutually  fully 
coherent  at  each  frequency.  These  are,  however,  exceptional  cases.  In  general, 
9r//Q((o)  will  be  frequency-dependent  and  Eq.  (3.4)  shows  that  the  field  spectrum  Sy{(xi) 
will  then  be  no  longer  proportional  to  the  source  spectrum  Sgfco).  Hence,  in  general, 
not  only  the  source  spectrum  but  also  the  correlation  between  the  two  sources  will 
determine  the  spectrum  of  the  emitted  light. 

Before  we  proceed  to  examine  the  effect  of  the  degree  of  spatial  correlation  on 
the  spectrum  of  the  radiation  we  simplify  our  notation  by  setting 

d2 

=  (3.7) 

We  win  refer  to  5v(CD)  as  the  reduced  field  spectrum.  For  the  sake  of  simplicity  we  will 
consider  source-correlations  that  are  characterized  by  a  real  degree  of  coherence. 
Equation  (3.4)  then  becomes 

rj^(F;<o)  =  5g(<u)[l-«'/ig(cy)]  .  (3.8) 

It  foUows  at  once  from  this  formula,  that  in  terms  of  sy  and  Sq, 


SyiCO)  ^ 

Sq(co) 


(3.9) 


From  the  inequality  (3.5)  and  from  Eq.  (3.9)  it  follows  that  only  those  reduced  field 
spectra  rv(w)  can  be  generated  for  which 
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s^{(0)<1Sq{0})  .  (3.10) 

5v<0))  and  Sg(a))are,  of  course,  necessarily  non-negative. 

Conversely,  when  the  inequality  (3.10)  is  satisfied  one  finds  at  once  from 
Eq.  (3.9)  that 


.  (3.11) 

Since  the  inequality  (3.1 1)  is  the  only  constraint  that  the  degree  of  spectral  coherence 
/ig(co)  must  satisfy  we  see  from  Eq.  (3.10) ,  that  any  reduced  field  spectrum  Sy{(ii) 
which  does  not  exceed  twice  the  magnitude  of  the  source  spectrum  5g(CD)  at  any 
frequency  to,  can,  in  principle,  be  generated  by  this  mechanism. 

We  now  demonstrate  how  source  correlations  can  give  rise  to  various  types  of 
spectral  changes. 


3.3  Change  in  spectral  linewidth 

Let  us  assume  first  that  the  source  spectrum  consists  of  a  single  spectral  line  of  a 
Lorentzian  profile,  viz. , 


S^(co)  = 


A) 

r^+(o)-Q)Qf  ’ 


(3.12) 


(coq,  Fq,  Aq  are  positive  constants  and  Fq  «  coq).  Suppose  that  we  wish  to  produce  a 
reduced  field  spectrum  that  consists  also  of  a  line  of  Lorentzian  profile  centered  on  the 
same  frequency  cOq,  but  is  of  a  different  width  and  of  different  strength,  say 
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S\/(^CO')  —  ~  H  ''  y 


(3.13) 


(  Fj,  Aj  arc  positive  constants  and  Fj  «  0)q  ).  On  substituting  ttom  Eqs.(3.12)  and 

(3.13)  into  the  inequality  (3.10)  we  find  that  we  must  have 


(3.14) 


where  [/()(«»)]  is  the  maximum  value,  in  the  range  0  <  0)<  »,  of  the  function 


r^+((o-o)Qr 


(3.15) 


Straightforward  calculation  shows  that  for  all  (positive)  frequencies  to. 


(ro/r,fs /<,(«)<! 

when 

r^>  Fq  , 

i</„«B)<(r,/r,f 

when 

F,<Fo  . 

(3.16) 


Using  these  inequalities  we  deduce  at  once  from  (3.14)  that  we  must  have 
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Relative  Frequency 

Figure  3.2:  An  example  of  line  broadening  by  source  correlations.  The  source  spectrum  (solid  line) 
and  the  reduced  field  spectrum  (dashed  line)  are  lines  of  Lorentzian  profiles  and  the  curves  are  normalized 
so  that  the  source  spectrum  has  the  value  unity  at  the  center  frequency.  The  relative  linewidths  used  are 
ro/(oo  =  o.i,  ri/(Do  =  0.15. 


<  2  when  r,>  Fr,  , 

(3.17) 

.^S2{r./r„f  when  r,<r„. 


In  the  first  case  (Fj  >  Fq)  the  emitted  (reduced)  spectral  line  is  broader  than  the 
spectral  line  of  the  source;  in  the  second  case  (Fj  <  Fq)  it  is  narrower. 

With  the  conditions  (3.17)  assumed  to  be  satisfied,  the  degree  of  spatial 
coherence  that  will  produce  the  reduced  field  spectrum  (3.13)  from  the  source  spectrum 
(3.12)  is  obtained  at  once  on  substituting  from  these  equations  into  the  formula  (3.9). 
One  then  finds  that  the  required  degree  of  spatial  coherence  is  given  by 
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Figure  3.3:  An  example  of  line  narrowing  by  source  cwielations.  The  source  spectrum  (solid  line) 
and  the  reduced  field  spectrum  (dashed  line)  are  lines  of  Lorentzian  profiles  and  the  curves  are  normalized 
so  that  the  source  spectrum  has  the  value  unity  at  the  center  frequency.  The  relative  linewidths  used  are 
To/mo  *  0-2.  r  i/tOQ  =  0.1. 


Aq  r^  +  ico-coQ^ 


(3.18) 


In  Figs.  3.2  and  3.3  wc  show  examples  of  line  broadening  and  line  narrowing 
with  Lorentzian  lineshapes.  The  spatial  correlation  that  produces  the  reduced  field 
spectrum  depicted  in  Fig.  3.3  is  plotted  as  a  function  of  frequency  in  Fig.  3.4. 

The  method  of  producing  Lorentzian  field  spectra  of  a  controlled  linewidths  can 
also  be  applied  to  Gaussian  lineshapes.^  The  source  spectrum  then  has  the  form 

Sg(cd)= /V  exp[-(G)-o)(,)^/25j]  .  (3.20) 


Spectral  modulation  with  scalar  sources 


Chapters 


79 


Relative  Frequency 

Figure  3.4:  The  degree  of  spatial  coherence  giving  rise  to  the  reduced  field  spectrum  of  a  narrower 
profile  shown  in  Fig.  3.3. 

(<%  Aq  axe  positive  constants  and  Sq  «  (Oq).  We  now  consider  the  possibility  of 
producing  a  reduced  field  spectrum  that  also  consists  of  a  line  of  Gaussian  profile, 
centered  at  the  same  frequency  (i\)  but  of  different  width  and  different  strength,  say 

5^,(0)  =  i4jexp[-(a)-cOo)^/25,^]  (3.21) 

(5i,i4i  are  positive  constants  and  «  cOq)-  O"  substituting  from  Eqs.(3.20)  and 
(3.21)  into  the  inequality  (3.10)  we  deduce  at  once  that  the  following  condition  must 
be  satisfied: 

,  (3.22) 
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where  [g(co)]num  is  the  maximum  value,  in  the  range  0  <  co<  «>  of  the  function 


g(tu)  =  exp[-(<a-n)o)^/i4]  , 


(3.23) 


where 


1  _  1 
A  ^ 


(3.24) 


Consider  first  the  case  when  5i  <  ^  (line  narrowing).  In  this  case  d  >  0  and 
evidendy  [g(<o)]m«x=  g(fJo)  =  1-  Hence  the  realizability  condition  (3.22)  becomes 


VA>^2. 


(3.25) 


On  the  other  hand  when  5i  >  ^,  d  becomes  negative  and  g{(0)  has  then  no  upper 
bound  in  the  range  0  <  oo<  oo .  Hence  a  broader  line  of  Gaussian  profile,  centered  at 
the  same  frequency  coj),  cannot  be  produced  by  this  mechanism.  However,  in  practice 
one  is  unlikely  to  be  interested  in  situations  where  the  spectra  Sq{(0)  and  Sy{(0)  have 
Gaussian  forms  for  all  frequencies.  If  one  requires  that  the  reduced  field  spectrum  has 
a  Gaussian  shap<*  only  over  a  finite  range  around  frequency  fflh*  say 


(OQ-a^co^Q)Q+P  , 


(3.26) 


where  a  and  p  are  positive  constants,  the  inequality  (3.22)  needs  only  be  satisfied 
when  the  maximum  of  g(a)  is  taken  over  the  restricted  range  (3.26).  Instead  of  the 
inequality  (3.22)  we  then  have  the  constraint 


Spectral  modulation  with  scalar  sources 


Chapters 


Relative  Frequency 

Figure  3.5:  An  example  of  line  narrowing  by  source  correlations.  The  source  spectrum  (solid  line) 
and  the  reduced  field  spectrum  (dashed  line)  are  lines  of  Gaussian  profiles  and  the  curves  ate  normalized  so 
that  the  source  spectrum  has  the  value  unity  at  the  center  frequency.  The  relative  linewidths  used  are 
Sq/Wo  ”  0-05,  Sj/(Dq  =  0.01. 

^exp(y2/2|4|)<2,  (3.27) 

A) 

where  y  is  the  largest  of  the  constants  a,  p. 

Returning  to  the  first  case  (A  <  Sq),  the  degree  of  spatial  coherence  needed  to 
achieve  this  modification  of  the  spectral  line  is  according  to  Eqs.(3.9),  (3.20)  and 
(3.21)  given  by 

MQi(0)  =  -^exp[-(<o  -  - 1  ,  (3.28) 

A) 
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Figure  3.6:  An  example  oi  line  broadening  by  source  correlations.  The  source  spectrum  (solid  line) 
and  the  reduced  field  spectrum  (dashed  line)  are  lines  of  Gaussian  profiles  and  the  curves  are  normalized  so 
that  the  source  spectrum  has  the  value  unity  at  the  center  frequency.  The  relative  linewidths  used  are 
~  Sj/{Oq  =  0.7. 


where  A  is  defined  by  Eq.(3.24).  In  the  second  case  (5i  >  6^)  the  degree  of  spatial 
coherence  is  given  by  Eq.(3.28)  only  for  frequencies  that  are  within  the  range  (3.26); 
for  frequencies  outside  this  range  the  degree  of  spatial  coherence  can  take  on  arbitrary 
values,  subject  to  the  constraint  expressed  by  Eq.  (3.11).  Line  narrowing  and  line 
broadening  of  Gaussian  profiles  are  illustrated  in  Figs.  3.5  and  3.6. 


3.4  Spectral  lineshifts 

We  now  consider  the  possibility  of  changing  not  only  the  width  of  the  line  but  also  its 
center  frequency,  from  (Oq  to  o)i  say.  Suppose  that  the  source  spectrum  5|2(^)  again 
the  single  spectral  line  (3.13)  of  Lorentzian  profile,  but  that  the  reduced  field  spectrum. 
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Figure  3.7:  Spectral  shift  and  tine  narrowing  with  Lorentzian  profiles.  The  source  spectrum  (solid 
curve)  has  a  relative  linewidth  rQ/(Oo  =  0.05  giving  rise  to  a  reduced  field  spectrum  (dashed  line)  of 
relative  linewidth  F j/Mq  =  0.01  centered  at  frequency  to  =  1  .ItOQ. 


whilst  also  a  line  of  Lorentzian  profile,  is  centered  at  a  different  frequency  (Oq,  i.e. 
that 


5v,(<o)  = 


(3.29) 


(o)|,  F],  Ai  are  positive  constants,  Fj  «  tOj).  On  substituting  from  Eqs.  (3.12)  and 
(3.29)  into  Eq.  (3.10)  we  find  that  the  following  inequality  must  now  be  satisfied: 


(3.30) 


Here  [fi (©)],„„  is  the  maximum  value  in  the  range  0  <  to  <  «>  of  the  function 
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Figure  3.8:  The  degree  of  correlation  that  gives  rise  to  the  firequency-diifted  reduced  field  spectrum 
shown  in  Fig.  3.7. 


/l(^0)  = 


rQ+jw-apf 

+ {to- co^)^ 


(3.31) 


Unlike  in  the  case  considered  above  (when  tU]  =  Wq)  an  explicit  expression  for  the 
maximum  value  of  this  function  cannot  readily  be  obtained  (see  Appendix  A  for  details 
of  the  numerical  method  used  to  determine  the  maximum  value  of  Aj).  However  it 
seems  that  with  suitable  choices  of  the  constants  that  specify  the  reduced  field  spectrum 
(3.29)  the  inequality  (3.30)  can  be  satisfied  for  all  non-negative  frequencies  ©.  The 
degree  of  spatial  coherence  which  gives  rise  to  the  reduced  field  spectrum  (3.29)  is  then 
obtained  on  substituting  from  Eqs.f3.29)  and  12^'  mto  Ea.(3.9V  Tf.c  resulting 
expression  is 
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/^(2(«a)  = 


.4^  r^+ico-co^f 


(3.32) 


In  Fig.  3.7  we  show  an  example  of  spectral  changes  in  which  both  the  linewidth  and 
the  center  frequency  are  modified  by  source  correlations  which  are  represented  by  the 
degree  of  spatial  coherence  in  Eq.  (3.32).  Figure  3.8  shows  the  degree  of  spatial 
coherence  that  gives  rise  to  these  spectral  changes. 


3.5  Modulation  of  a  single  spectral  line 

As  a  last  example  we  consider  the  possibility  of  generating  from  a  source  spectrum  that 
consists  of  a  single  line  of  Lorentzian  profile  a  field  spectrum  that  consists  of  several 
lines  of  Lorentzian  profile.  More  specifically,  with  two  sources  that  have  identical 
spectra  given  by  Eq.  (3.12)  we  wish  to  generate  a  field  whose  reduced  spectrum  has 
the  form 


N 


Syi(0)-  y  ^ 


(3.33) 


where  N,  and  r„  are  positive  constants  and  r„  «  (1  <  /i  ^  A^) .  For  this 

to  be  possible  the  following  condition  obtained  on  substituting  Eqs.  (3.33)  and  (3.12) 
into  Eq.  (3.10)  must  be  satisfied: 
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Figure  3.9:  Generation  of  three  spectral  lines  from  a  single  spectral  line,  all  of  Lorentzian  profiles. 
The  source  spectrum  consists  of  a  single  line  whose  relative  width  is  TQ/tao  =  02S  and  the  three  spectral 
lines  are  all  of  relative  width  F/tOo  »  O.OS. 


1 


max 


<2  . 


(3.34) 


Here  [fn(o))]max  is  the  maximum  value,  in  the  range  0  ^  <  «>o,  of  the  function 


(3.35) 


Assuming  that  the  constraint  (3.34)  is  satisfied,  the  degree  of  spatial  coherence  needed 
for  generating  the  reduced  field  spectrum  (3.33)  is  obtained  at  once  on  substituting 
from  Eqs.  (3.33)  and  (3.12)  into  Eq.  (3.9).  One  then  finds  that 
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Figure  3.10:  The  spatial  correlation  for  generation  of  three  spectral  lines  from  a  single  spectral  line. 


,  .  I  ^  .  r'n+io)- 


r^+jco-coQf 


-1 


(3.36) 


An  example  of  generation  of  a  field  spectrum  consisting  of  three  lines  from  a 
source  spectrum  consisting  of  a  single  line  is  illustrated  in  Fig  3.9.  Figure  3.10  shows 
the  degree  of  correlation  which  gives  rise  to  this  spectrum. 


3.6  The  spectrum  produced  by  an  array  of  partially  coherent 
sources 

So  far  we  have  considered  a  very  simple  radiating  system  consisting  of  two  small 
sources  with  identical  spectra  and  we  showed  that  by  appropriately  correlating  them. 
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the  spectrum  of  the  tmitt-d  radiation  can  take  on  many  different  forms.  With  systems 
consisting  of  a  larger  number  of  radiating  sources  one  can,  of  course,  produce  more 
diverse  spectral  changes  than  those  considered  above.  A  simple  extension  of  the 
system  which  we  considered  is  an  array  of  2/V  equally  spaced  small  fluctuating 
sources^  (see  Fig.  3.11).  The  field  produced  by  the  array  is  given  by  [cf.  Eq.  (3.2)) 

N  ikR. 

V(F;co)=^(2.(a))^.  (3.37) 


Using  Eq.  (3.37),  the  spectrum  of  the  field  at  the  point  P  is  readily  found  to  be 


Sy{P‘,(0)  —  Sq{(0) 


N 


1 

i--N  i*j 


R,Rj 


(3.38) 
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Equation  (3.38)  shows  some  interesting  new  features.  Even  for  observation  points  on 
the  axis,  the  exponential  factor,  does  not,  in  general,  reduce  to  a  constant  which  is 
independent  of  frequency.^  As  a  result,  it  follows  that  the  spectnim  of  the  field  will  no 
longer  be  proportional  to  the  spectrum  of  the  sources  when  the  fluctuations  are  fully 
correlated,  /iij(ty)  =  1.  Specifically,  the  presence  of  the  wavenumber  k  in  the 
exponential  causes  some  modulation  of  the  spectrum  that  is  proportional  to  the 
magnitude  of  the  degree  of  corrclation.^  On  the  other  hand,  when  the  source 
fluctuations  are  uncorrelated,  =  0  when  i  ^  j  and  we  then  have 

2N 

5v,(P;tD)  =  Sg(0))]^/?j  ^  ,  (3.39) 

i=l 

The  spectrum  of  the  field  is  now  seen  to  be  proportional  to  the  spectrum  of  the  source 
for  all  points  of  observation. 


3.7  General  remarks  on  spectral  modulation 

In  thi  chapter  we  presented  simple  examples  of  partially  coherent  radiating  systems  in 
which  the  field  spectrum  is  modified  by  the  degree  of  correlation  of  the  source 
fluctuations.  As  is  evident  from  Eqs.  (3.3)  and  (3.38),  this  mechanism  does  not 
produce  any  new  ficquencies  but  just  increases  ot  decreases  their  relative  contributions 
to  the  observed  spectrum.  This  fact  implies  that  in  order  to  generate  field  spectrum  of  a 
particular  shape,  the  source  must  possess  all  the  required  spectral  components.  For  this 
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reason,  it  is  convenient  to  choose  broad  t  .nd  sources  tor  theoretical  and  experimental 
illustration  of  the  spectral  effects  which  we  discussed  in  this  thesis. 

The  basic  question  that  we  have  not  addressed  in  this  chapter  is  how  to  produce 
the  prescribed  spaMal  cc^relations.  Several  methods  for  generating  and  modifying 
source  correlations,  at  least  for  secondary  sources,  have  been  developed  in  recent 
years.  They  include  the  use  of  scattering  bv  liquid  crystals  under  the  influence  of  an 
external  D.C.  field^,  the  use  of  rotating  ground  glass  plates*  and  of  holographic 
filters^,  interaction  of  light  with  ultrason’c  wavesi^by  imaging  -^Jid  lensless  feedback 
systems  and  by  t!ie  use  of  achromatic  Fourier  transform  lenses. n 

The  most  notable  method  and  experimental  demonstration  of  spectral 
modulation  by  control  of  source  correlation  was  reported  by  G.  Indebetouwi*.  His 
technique  makes  use  of  spatial  masks  which  produce  a  secondary  source  of  prescribed 
frequency  dependent  degree  of  correlation. 
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spectrum  away  from  line  center. 
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Appendix  A:  The  method  used  in  determining  the 

maximum  value  of  Aj 

Starting  from  Eqs.  (3.30)  and  (3.31)  we  have  the  constraint 


- * - - - TT  S  /  . 

Aq  ^^+(co-o)^r 


First  we  set 


A.  =  <A  2) 

to  normalize  the  source  spectrum  to  the  value  of  unity  at  co  =  (Dq-  Next  we  set  the 
frequency  ©  =  ©j  in  Eq.  (A  1).  This  choice  is  made  only  for  computation  purposes  and 
should  not  be  mistaken  as  a  statement  that  the  function  fi(©)  in  Eq.  (3.30)  is  maximized  at 
©  =  ©1.  It  follows  from  Eqs.  (A  1)  and  (A  2)  that  the  upper  bound  on  Aj  in  this  case  is 
given  by 


[A]  = 

L  ^Jmax 


r^+(o)i-Q)of 


(A3) 


This  method  for  determining  the  value  for  A]  was  used  with  numerous  choices  of 
the  constants  Fq,  Tj,  ©q,  ©i  and  the  resulting  degree  of  correlation  has  alwas  satisfied  the 
constraint  (3.1 1). 
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Introduction 

The  effects  of  source  correlations  on  the  spectrum  of  emitted  radiation  have  been  treated 
so  far  only  in  the  framework  of  scalar  theory.  In  the  previous  chapter  we  considered  a 
single  system  consisting  of  two  scalar  sources  and  we  examined  some  of  the  spectral 
effects  that  may  be  produced  on  the  axis  by  controlling  the  correlations  between  the  two 
sources.  In  this  chapter  we  consider  a  system  of  two  partially  correlated  linear  dipoles. 
First,  we  derive  expressions  for  the  observed  spectrum  at  an  arbitrary  point  in  the  far 
zone.  We  then  examine  how  the  degree  of  correlation  of  the  dipole  polarization 
fluctuations  affects  the  spectrum  of  the  field  observed  in  particular  directions.  We  also 
consider  the  effects  of  the  degree  of  correlation  on  the  angular  distribution  of  the 
radiant  intensity  for  fixed  frequencies.  In  Secs.  4.4  and  4.5  we  present  expressions  for 
the  total  radiated  power  and  for  the  directivity  of  the  system  as  a  function  of  the  degree 
of  correlation  and  the  spatial  separation  between  the  two  dipoles.  We  also  compare  our 
results  for  the  partially  correlated  dipoles  with  results  in  the  well  known  limiting  cases 
of  fully  correlated  and  uncorrelated  dipoles.^  This  comparison  gives  a  valuable 
measure  for  the  range  of  possible  spatial  and  spectral  modulation  effects  which  can  be 
produced  by  controlling  source  correlations. 


4.2  Far  zone  spectrum  of  partially  correlated  dipoles 

Consider  two  linear  dipoles,  situated  at  points  ±  yQ  and  vibrating  in  the  z  direction  as 
shown  in  Fig.  4.1. 

Let 
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P^{r,t)  =  p^{t)5ir-y^)z, 

P2ir.t)  =  P2{t)5{r-\-y^)z, 


(4.1) 


be  the  electric  polarization  vectors  of  the  two  dipoles.  Here  p^{t),  (j  =  1,2),  specify 
the  polarization  fluctuations  of  the  dipoles  as  functions  of  time  and  ^  and  i  are  unit 
vectors  in  the  positive  y  and  z  directions  respectively.  We  assume  that  pjj)  are  random 
functions  of  time,  characterized  by  stationary  ensembles. 

In  the  space-frequency  representation,  the  electric  Hertz  potential  of  the  field 
produced  by  the  dipoles  is  given  by 


n^(r;o))  =  z 


ikJt, 


UcR^ 


Pl(C0)—  +  P2{(0)—- 


(4.2) 


where  -  yoy|,  ^  =  |r  +  >(,>1 ,  and^ 


z 


Figure  4.1:  Illustrating  the  configuration  and  the  notation.  The  two  dipoles 
are  marked  by  the  heavy  vertical  arrows  are  separated  by  a  distance  2yQ  bom  each 
other. 


Spectral  effects  with  electromagnetic  sources 


Chapter  4 


97 


Pjico)=  \pjit)e‘^dt  0  =  1,2).  (4.3) 

— oo 

In  order  to  calculate  the  spectral  intensity  in  the  far  zone,  it  is  sufficient  to 
evaluate  the  magnetic  field  only;  it  is  given  by  the  expression: 

B{r,0})  =  -ijtV  X  n^ir;a)).  (4.4) 

Since  n^(r\co)  is  a  vector  along  the  z-direction,  it  follows  that 


Vxn^(r;(o)=  4-x--^y  n{r,co). 

[dy  ox  j 


On  substituting  Eq.  (4.2)  in  Eq.  (4.5)  and  calculating  the  partial  derivatives,  we  obtain 
the  formula 


Vxi^(r,<B)  =  pj(co)-—  ik-  — 


L  Z-Zo-  Ar. 


rJ[  R,  ^  r/ 


^  1  Z+Zo.  ^ 


For  field  points  r  =  ru  in  the  far  zone  kR^  »  1  (j  =  1,  2)  and  we  have 
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Figure  4.2:  Gr^hical  rqnesentadon  of  the  approximation  fcM’/?,  made  in  Eq.  (4.7). 


/2l~r-yo  («•>), 

R^-^r  +  y^iu-y), 

where  M  is  a  unit  vector  in  the  direction  of  observation  (see  Fig.  4.2). 

Using  Eqs.  (4.6)  and  (4.7),  we  may  express  the  far  zone  5-field  in  the  form 

B{ru\w) — - sin  6 

r 

[pj  (Q>)exp(-i7:yo«  •  y)  -i-  p^{G))txp[ikyQU  •  y)]  (4.8) 

as  kr-^oo. 

Here,  we  have  used  the  spherical  polar  coordinates  (r,  6, 4>)  with  the  polar  axis  along 
the  z-direction  and  with  0  =  sin  6(y  cos^  - xsin  ^). 

The  radiant  intensity  J{u\(o) ,  i.e.  the  power  per  unit  frequency  at  frequency  on, 
per  unit  solid  angle  around  the  direction  specified  by  the  unit  vector  u,  is  given  by 
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J{u\(o)=  lim  P'u 


(s^'"Hru;o))y 


(4.9) 


where  S^‘‘\ru;Q))  is  the  Poynting  vector  at  points  in  the  far  zone,  and  the  angular 
brackets  denote  the  ensemble  average.  In  terms  of  the  B-field,  we  have 


J(u;(i))=  Um  x  iij  xJ8^“\r«;to)  ^  (4.10) 


where  SR  denotes  the  real  part  and  c  is  the  vacuum  speed  of  light  Using  the  vector 
identity 


ll■(J?XIl)*  xB  =  |i?Xttp,  (4.11) 

we  simplify  Eq.  (4.10)  for  the  radiant  intensity  to  the  form 

/(«;(»)=  lim  r^-^(|J?(m;a))|^y  (4.12) 

kr— »«>  o7l  '  ' 


Let  us  now  assume  that  both  dipoles  have  the  same  spectrum,  i.e.,  that 

(lPi«u)f )  =  =  S^ico).  (4.13) 

We  also  introduce  the  complex  degree  of  spatial  coherence  at  frequency  to,  which 
characterizes  the  correlation  between  the  two  dipoles,  by  the  formula 
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^lp{(0)^{p^{(0)p^{O)))IS^{(0).  (4.14) 

On  substituting  from  Eqs.  (4.8),  (4.13)  and  (4.14)  in  Eq.  (4.12),  we  obtain  the 
following  expression  for  the  radiant  intensity  produced  by  the  two  partially  correlated 
dipoles: 


,4 

J(u;(o)  =  5p(6))sin^  ^jl  +  9l^/i^(fi))exp(2//tyQ  sin  0  sin  ^)j|.  (4.15) 


If  we  express  the  degree  of  correlation  in  the  form 


p„((o)  =  \^p(wY 


Vt0») 


(4.16) 


where  2yf{ci))  corresponds  to  the  effective  steering  angle^  ,  we  find  that  the  radiant 


intensity  of  the  two  partially  correlated  linear  dipoles  is  given  by  the  formula 


ck^ 

Jiu;  (0)  =  ^^5^(m)sin^  ©jl  +  |jU^(a)^cos[2ily(,  sin  0  sin  0  +  2  V'^(Q))]|. 


(4.17) 


4.3  The  elTects  of  spatial  correlation  on  the  spectrum  and 
the  angular  distribution  of  the  radiant  intensity 

We  now  consider  some  special  cases  which  will  help  to  illustrate  the  significance  of  Eq. 

(4.17) .  First  let  us  consider  two  uncorrelated  dipoles.  In  this  case  Pp(.co)  =  0,  and  Eq. 

(4.17)  reduces  to 
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Ck*  2 

= - 5  (G))sin^  0  . 

4n  ^ 


(4.18) 


As  one  may  expect,  the  same  spectrum  is  obtained  if  the  radiation  originated  from  a 
single  dipole  located  at  the  origin  and  having  polarization  fluctuations  equal  to  the  sum 
of  the  polarizations  of  the  two  dipoles. 

Similarly,  when  the  two  dipoles  are  fully  correlated,  i.e.  when 

the  formula  (4.17)  for  the  radiant  intensity  reduces  to 

=  ~Sp(ft))sin^  0cos^[*yo  sin  Osin  (j>  +  V^(G))].  (4.20) 


(4.19) 


We  observe  that  for  every  frequency  to,  the  phase  angle  V'(t>))  and  the  separation 
constant  I^q  completely  determine  the  angular  distribution  of  the  radiated  power. 

Returning  to  the  general  case  when  the  two  dipoles  are  partially  correlated,  we 
note  that  when  the  point  of  observation  is  on  the  x-axis  (6=  a/2,  <f>  =  0),  the  radiant 
intensity  according  to  Eq.  (4.17)  is  given  by 

Jix-,0))  =  (4.21) 

The  result  for  this  special  case  is  in  the  same  form  as  the  corresponding  expression  for 
the  radiant  intensity  from  two  small  partially  correlated  scalar  sources  [cf.  Eq.  (3.6)] 
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Sy{P-,0})  =  ^SQ{(o\\  +  ^.^Q{(o)\  (4.22) 

where  P  is  any  point  on  the  perpendicular  bisectw  of  the  line  adjoining  the  two  sources. 
Aside  from  the  2//?^  factor,  the  only  significant  difference  is  in  the  factor  appearing 
in  Eq.  (4.21),  and  as  a  result,  the  product  k*S^((0)  is  shifted  to  higher  frequencies 

relative  to 

The  radiation  pattern  generated  by  the  partially  correlated  dipoles  differs  from 
the  two  limiting  cases  of  fully-correlated  and  of  completely  uncorrelated  dipoles  in 
several  ways.  If  we  denote  the  direction  for  which  the  radiant  intensity  is  a  maximum 
by  11^,  it  can  readily  be  shown  that  the  maximum  possible  radiant  intensity 

[y(i/^;ci))]^^  produced  by  the  two  dipoles  is  obtained  when  they  are  fully  correlated, 
and  is  given  by 


2jc 


(4.23) 


On  the  other  hand  when  the  dipoles  are  only  partially  correlated,  the  maximum  radiant 
intensity  is  smaller  by  a  factor  ^|l  +  |/i^(a))|J. 

In  the  two  limiting  cases  of  fully-correlated  dipoles  and  uncorrelated  dipoles, 

A 

the  nulls  of  the  radiant  intensity  distribution  are  determined  by  the  factor  sin  0  in 
Eq.  (4.18),  and  by  the  factor  sin^  0cos^[ilyQsin0sin0+  V'^(O))]  in  Eq.  (4.20).  By 
contrast,  it  follows  from  Eq.  (4.17)  that  when  the  dipoles  arc  partially  correlated, 
0  <  <  1,  there  are  no  nulls  outside  the  plane  0  =  0.  This  fact  is  significant  in 

connection  with  the  theorems  regarding  the  approximation  of  desired  radiation  patterns 
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by  arrays  of  radiators.^  Specifically,  in  cases  where  the  directivity  of  antennas  is 
important,  effort  is  made  to  reduce  the  amount  of  sidelobes.  Typical  antennas  in  which 
all  elements  are  radiating  coherently  have  several  nulls  in  their  radiation  pattern  and 
hence  there  is  a  latitude  of  design  options  allowing  for  the  shifting  of  these  nulls.  On 
the  other  hand,  the  system  considered  here  has  no  nulls  in  the  radiation  pattern  outside 
the  plane  0  =  0  in  which  the  radiation  vanishes  identically.  As  a  result,  no 
rearrangement  of  the  elements  can  produce  a  null  in  a  particular  direction  as  long  as  the 
degree  of  correlation  is  less  than  unity. 

We  now  illustrate  the  angular  distribution  of  the  radiant  intensity  produced  by 

the  two  partially  correlated  dipoles.  Let  the  dipole  spectrum  be  a  Lorentzian  line  of 
width  centered  at  frequency  co^,  i.e.. 


5/0))  = 


1 

T+OoT^’m^TT^' 


(4.24) 


We  choose  a  real-valued  degree  of  correlation  in  the  form 
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Relative  Frequency 

Figure  4.3:  The  degree  of  correlation  given  in  Eq.  (42S).  The  constants  used  arc  (0,/(Oo  =  0.8. 
tOj/tOfl  =  1.2,  8,  =  Sj  =  0.05  and  A,  =  Aj  =  1.98. 

Hp((0)  =  \  exp[-(<y  -  4  exp[-(©  -  (4.25) 

Here  Cdj,  (O^,  Sy  and  arc  positive  constants,  and  we  select  and  so  that  the 
constraint 

|/i/a»|^l.  (4.26) 

is  satisfied  throughout  the  frequency  range  of  interest  We  note  that  by  choosing  a  real- 
valued  degree  of  correlation  the  steering  angle  is  set  to  zero  (i.e.,  2v/(a))  =  0)  in  the 
xollowing  examples. 
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Figure  4.4:  Angular  distribution  of  the  radiant  intensity  at  frequency 
(U  =  (Og  for  fully  correlated  dipoles  (solid  line)  and  partially  correlated 
dipoles  (dashed  line).  The  separation  constant  is  ityg  =  1  and  the  degree  of 
correlation  is  given  by  Eq.  (4.25). 


Figures  4.4  and  4.5  show  the  angular  distributions  of  the  radiant  intensity  in  the 
X,  y-plane  (9  =  rr/2)  for  two  fully  correlated  dipoles  and  foi  partially  correlated 
dipoles  whose  degree  of  correlation  is  given  by  Eq.  (4.25).  In  Fig.  4.4  the  angular 


0»  ISC' 


Figure  4.5:  The  angular  distribution  of  the  radiant  intensity  for  fully  correlated  dipoles  (solid  line) 

and  partially  correlated  d’poles  (dashed  line)  and  separation  constant  lcyo=\.  The  relative  frequencies 
are  (o/Wg  =  1.05  [(a)]  and  0)/(i)g  =1.1  [(b)]. 
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distribution  is  calculated  for  the  center  frequency  CD  =  of  the  polarization  spectrum, 
while  Fig.  4.5  shows  the  angular  distributions  at  different  frequencies. 

We  note  that  in  spite  of  the  fact  that  the  steering  angle  was  fixed,  r»Kco)  =  0], 
there  is  a  directional  shift  in  the  angular  distribution  of  the  radiant  intensity  produced  by 
partially  correlated  dipoles. 

When  the  separation  between  the  dipoles  increases  the  lobe  structure  becomes 
more  complicated.  In  Figs.  4.6  and  4.7  we  show  the  angular  distributions  of  the 
radiant  intensity  for  dipole  separation  =  3.  Comparing  the  angular  distributions 
for  dipole  separation  =  1  (Fig.  4.4,  4.5)  and  k^Q  =  3  (Fig.  4.6,  4.7),  we  note 
that  at  each  frequency  the  number  of  lobes  with  our  two  choices  of  correlations  is 
identical  although  their  angular  distributions  are  somewhat  different.  In  fact,  except 
when  the  two  dipoles  arc  uncorrelated,  the  number  of  lobes  of  the  radiant  intensity,  n, 
is  independent  of  the  degree  of  correlation  and  it  is  given  by 


n  = 


(4.27) 


90“ 


Figure  4.6:  The  angular  disuibution  of  the  radiant  intensity  for  fully 
correlated  dipoles  (solid  line)  and  partially  correlated  dipoles  (dashed  line)  and 
separation  constant  kyg  =  3.  The  relative  frequency  is  (o/Oq  =  1 . 
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180“ 


90“ 


0“ 


(a) 


270“ 


0“ 


Figure  4.7;  TTie  angular  distribution  of  the  radiant  intensity  fcM'  fully  correlated  dipoles  (solid  line) 
and  partially  correlated  dipoles  (dashed  line)  and  separation  constant  ky^  =  3.  The  relative  frequencies 
are  to/coo  =  1.05  [(a)]  and  (o/tOo  =  1.1  [(b)]. 


Here  the  square  brackets  denote  the  largest  integer  which  is  smaller  or  equal  to  the 
quantity  in  the  brackets. 

The  two  limiting  cases  of  fully  correlated  and  uncoirelated  dipoles  constitute  the 
limits  of  possible  modification  of  spectra.  In  Fig.  4.8  we  illustrate  the  range  of 
modulation  that  can  be  achieved  at  every  frequency  (o,  by  the  variation  of  the  magnitude 
of  the  degree  of  correlation  and  the  phase  2v^a>).  The  figtue  shows  four 

concentric  circles,  (/)  to  (iv),  representing  the  relative  ranges  of  the  angular  distributions 
of  the  radiant  intensity  in  the  x,y-plane.  The  external  circle  (0  corresponds  to  the  limit 
of  fully-correlated  dipoles.  As  one  may  observe  from  Eq.  (4.20),  the  radiant  intensity 
of  two  correlated  dipoles  may  attain  any  value  inside  the  circle  of  radius  ck^SJ,0))lliL 
The  locaticMi  of  the  maximum  of  the  radiant  intensity  is  determined  by  the  choice  of  the 
steering  angle  2y/ia}).  Similarly  it  follows  from  Eq.  (4.19)  that  when  the  dipoles  are 
unconelated,  the  radiant  intensity  is  consumned  to  the  uniform  value  on  a  circle  (Hi)  of 
radius  ck*S^i(o)/47[.  When  the  dipoles  are  partially  correlated  we  see  from  Eq.  (4.17) 


Spectral  effects  with  electromagnetic  sources 


Chapter  4 


108 
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Figure  4.8:  The  range  of  possible  values  of  the  angular  distribution  of  the 
radiant  intensity.  The  shaded  area  indicates  the  region  in  which  the  maxima  and 
minima  of  the  radiant  intensity  for  partially  correlated  dipoles  are  found. 


that  the  radiant  intensity  can  have  any  value  in  an  annular  domain  bounded  by  circles 
(I'O  and  (iv).  It  is  also  apparent  from  this  representation  that  unless  the  two  dipoles  are 
fully  correlated  there  arc  no  nulls  of  the  radiant  intensity  outside  the  plane  ©  =  0. 

Figure  4.8  gives  a  symbolic  representation  of  the  directivity  of  the  system 
consisting  of  two  partially  correlated  linear  dipoles.  A  detailed  discussion  of  the 
directivity  is  presented  in  Sec.  4.5. 


4.4  The  total  emitted  power 

The  total  power  Piw)  radiated  by  the  system  at  frequency  a>  is  given  by  the  expression 


P{o))=  jj{u;0))d£2, 

(471) 


(4.28) 
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where  the  integration  extends  over  the  whole  4ji  solid  angle.  On  substituting  from  Eq. 
(4.17)  in  Eq.  (4.28)  we  obtain  the  expression 


Pi(0)  =  ^^Spict))  jsin^0|l  +  9t^/i^(ci))exp(j2iyQsin0sin0)j|di2.  (4.29) 

4tt  , . _ . 


(4k) 


After  performing  the  integration  in  Eq.  (4.29)  [see  Appendix  B],  we  find  that 

Pi(0)  =  jck^Spico){l  +  /(2^o)9l[/ip(©)]}.  (4.30) 


where 


Figure  4.9:  The  behavior  of  the  function  fiz)  in  Eq.  (4.31)  as  a  function  of  the  parameter 
z  =  Iky^. 
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and  JqU)  and  j\  (z)  are  spherical  Bessel  functions. 

The  Tnaximum  of  the  function /(r)  can  be  shown  to  occur  when  z  =  0  where  it 
has  tne  value  f{0)  =  1  (see  Fig.  4.9).  It  follows  that  the  upper  bound  for  the  total 
radiated  power  is 


(4.32) 


In  the  limiting  case  when  the  two  dipoles  are  uncorrelated,  the  total  radiated 
power  is  given  by 


On  comparing  the  total  power  of  two  partially  correlated  dipoles  and  uncorrelated 
dipoles  [f*((‘‘^)]uncorr  ^9*’  (4.30)  and  (4.33)  that  the  ratio  between  the 

power  emitted  by  two  partially  correlated  dipoles  and  two  uncorrelated  dipoles  is  given 
by 


P(CD) 


=  I  +  /<2^(i)*[m,«»)]- 


(4.34) 


Since  the  function /(z)  in  Eq.  (4.30)  decreases  rapidly  with  increasing  z,  it  is  evident 
that  the  total  power  emitted  by  two  dipoles  that  are  separated  by  a  distance  that  is  much 
larger  than  a  wavelength,  is  equal  to  the  total  radiated  power  from  two  uncorrelated 
sources.  In  other  words,  when  the  separation  between  the  two  dipoles  is  large 
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compared  to  the  wavelength  of  the  radiation,  no  correlation  effects  could  be  deduced 
from  measurements  of  the  total  emitted  power. 


4.5  The  directivity  of  two  partially  correlated  dipoles 

The  directivity  D(u;  o))  of  a  radiating  system  is  defined  by  the  ratio® 


D{u;  CO)  =  AtU  (u;  coyPiw) , 


(4.34) 


where  u  is  the  direction  of  observation.  We  are  particularly  interested  in  the  maximum 
directivity  D(a))  s  co)  that  occurs  in  a  particular  direction  u^.  It  follows  from 
Eqs.  (4.23)  and  (4.30)  that  the  directivity  of  the  two  partially  correlated  dipoles  is  given 


3  l+p_(a)J 

D(a))  =  4 - .  . 

2l  +  fi2kyQ)<3^\i^i(o)] 


(4.35) 


When  the  two  dipoles  are  uncorrelated,  |/i^(£())j  =  0  and  the  directivity  reduces  to  the 
value  D(o})  =  3f2,  as  expected.  Similarly,  when  the  dipoles  are  fully-correlated  the 
directivity  is  given  by 


Di(o)  = 


l  +  /(2Ayo)cos[2v(to)]  ’ 


(4.36) 


where  2y/{co)  denotes  again  the  phase  of  /ipfcu). 


Spectral  effects  with  elecbomagnetic  sources 


Chapter  4 


112 


Equations  (4.35)  and  (4.36)  also  indicate  that  when  the  two  dipoles  are 
separated  by  a  distance  that  is  much  larger  than  a  wavelength,  the  directivity  of  partially 
correlated  dipoles  approaches  the  value 


lim 

b'o-*' 


Z)(0))  =  3[l  +  |/i/o»|]/2. 


(4.37) 


while  the  directivity  of  fully  correlated  dipoles  approaches  the  value 


lim  [D(n))]^^  =3. 


(4.38) 


We  also  note  that  if  we  use  Eq.  (4.16)  in  Eq.  (4.35)  and  differentiate  with  respect  to  the 
magnitude  of  the  degree  of  correlation  m  =  j;i^(r»)j,  we  find  that 


<iP((0)^  3{l  +  /(2^o)cos[2v(a))]] 
^  2[1  +  /7i/'(2^q)cos[2v(©)] 


>  0, 


(4.38) 


which  implies  that  for  a  fixed  value  of  the  quantities  f{2kyQ)co^2yfiQ))],  the 
directivity  is  an  increasing  function  of  the  magnitude  of  the  degree  of  ccnrelation.  As  a 
result  it  is  clear  that  the  maximum  directivity  is  obtained  for  fully  correlated  sources. 


4.6  Summary 

In  this  chapter  we  considered  the  effect  of  spatial  correlations  on  the  spectrum  of  the 
field  produced  by  partially  correlated  linear  dipoles.  We  showed  the  changes  in  the 
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angular  distribution  of  the  radiant  intensity  and  the  changes  in  the  spectrum  in  a  Hxed 
direction  as  a  function  of  correlation  of  the  dipoles  polarizaticm. 

In  examining  the  total  emitted  power  and  the  directivity  of  the  system  we 
showed  that  tht  maximum  radiated  power  and  the  maximum  directivity  are  obtained 
when  the  dipoles  are  fully  correlated. 
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Appendix  B:  Evaluation  of  the  integral  in  Eq. 

the  total  emitted  power 


Wc  begin  with  Eq.  (4.29),  viz., 

P(o))  =  ^s/ci»[;,+/j], 


where 

I^=  Jsin^edQ, 

(4ii) 


I2  =  Jsin^09i^jip(<i))exp(i2ilyQsin6sin<j>)jdQ. 
(4k) 


The  value  of  the  first  integral  is  readily  found: 


2k  k 

/,=  -cos^  6)sin6d0 

0  0 

=  2nj(l-x^)dx  =  — . 

«  ^ 


The  second  integral  is  given  by 


2k  k 

/j  =  J  d<t>Jsin^0exp(i2i^Qsin0sin<|>)d0d(|>. 
0  0 


In  evaluating  the  second  term  wc  first  perform  the  integral^  over  (ft,  giving 


(4.29)  for 


(Bl) 


(B2) 


(B3) 


(B4) 


(B5) 
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iqi 

/j  =  4k  J sin^ eyoClJfcyo sin0)de.  (B  6) 

0 


Using  the  identity 


sin^  0  =  sin  6  -  sin  0cos^  0, 


(B7) 


Eq.  (B  6)  can  be  expressed  as  a  sum  of  two  integrals 
*/2  %(! 

/2  =  4k  J  sin67Q(2iiyQsin0)d0  +  4K  Jsin0cos^0/o2^o*“®)^- 
0  0 


The  first  integral  in  Eq.  (B  8)  has  the  valued 


m 

4k  J sin  6JQ{2kyQ  sin  6)dd  ~  An  jQilky^),  (B  9) 

0 


and  the  second  integral  in  Eq.  (B  8)  has  the  value 


it/2 

4k  j  sin  0cos^  sin  6)dd  =  4k 

0 


y'lU^o) 

2*yo 


(BIO) 


Using  Eqs.  (B  4),  (B  9)  and  (B  10)  in  Eq.  (B  1)  we  obtain  the  following 
expression  for  the  total  emitted  power 
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P{(0)  =  |c/:%(0)){l  +  fi2kyQ)^n^i(0)^, 


(BID 


where 


/(D  =  |[yo(2)-;,(z)/z], 


(B12) 


and  JqU)  and  j\iz)  are  spherical  Bessel  functions. 
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Introduction 

The  theory  presented  in  the  previous  chapters  used  model  sources  for  the  representation 
of  die  cross- spectral  density,  characterized  by  a  spectral  density  function  and  a  spatial 
coirelation  function.  In  this  chapter  we  investigate  spectral  effects  which  are  produced 
by  source  eomelations  using  coherent  mode  representation  This  approach  emphasizes 
the  role  of  spatial  correlations  in  generating  spectral  changes  on  propagation  of  partially 
coherent  light,  arid  it  clarifies  the  distinction  between  diffraction  and  correlation  effects. 

We  will  consider  the  class  of  statistically  stationary,  secondary,  partially 
coherent,  planar  sources,  whose  cross-spectral  density  function  consist  of  a  finite  sum 
of  Hermite-Gaussian  modes.  We  evaluate  the  spectral  changes  which  occur  on 
propagation  by  decomposing  the  source  into  its  coherent  modes  and  evaluating  the 
contribution  of  each  mode  to  the  observed  spectrum.  We  determine  the  field  spectra  for 
observation  points  both  in  the  far  and  in  the  near  zone.  A  comparison  of  the 
calculations  shows  that,  at  least  in  the  cases  considered  in  this  chapter,  most  of  the 
spectral  changes  seen  in  the  far  zone  arc  already  present  in  the  spectrum  that  would  be 
ob.'^erved  at  very  short  distances  from  the  source  plane.  A  detailed  discussion  of  the 
development  of  spectral  changes  with  propagation  distance  from  '.’le  source  plane  is 
given  in  Chapter  6. 

Since  we  are  considering  sources  consisting  of  Hermite-Gaussian  modes  and 
because  such  modes  also  represent  the  transverse  modes  of  certain  laser  cavities^,  our 
analysis  provides  an  indication  of  the  type  of  spectral  effects  might  be  expected  to  arise 
in  fields  produced  by  some  multi-mode  lasers. 
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5.2  Coherent-mode  representation  of  partially  coherent 
sources  and  fields 

The  cross-spectral  density  function  of  a  statistically  stationary,  planar,  secondary 
source  of  any  state  of  coherence,  occupying  a  finite  domain  D  may  be  expressed  as  a 
sum  of  coherent  modes  viz  3, 

=  SX,.„((0)4„(r,;<I.)x,,,(rj;oo).  (5.1) 


Here  Tj  and  are  position  vectors  of  points  in  the  source  domain  D,  ^ 

eigenvalues  and  ^  eigenfunctions  of  the  Fredholm  integral  equation 

A  =  (5.2) 

D 

Here  the  eigenvalues 


{alln,m).  (5.3) 

The  eigenfunctions  (modes)  of  *1*®  source  are  taken  to  be  orthonormal, 

i.e.. 


=  (5.4) 

D 

where  5  is  the  Kronecker  symbol. 

The  spectrum  of  the  field  at  a  typical  observation  point  r  is  given  by 
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S(r;o))  =  SX^«(co)|v^^(^;(l))|^  (5.5) 


where  V„^(r,to)  is  given  by 

=  J  X».«(''’;t^)G(r,r';a))dV,  (5.6) 

D 

and  G(r,  r';a))  is  the  Green’s  function,  which  characterizes  propagation  from  the 
source  point  r'  to  the  field  point  r. 

It  is  evident  from  Eq.  (5.5)  that  each  term^  contributes  to  the  total  spectrum 
independently  of  all  the  others,  i.e.  there  are  no  cross-terms.  Moreover,  since  each 
term  under  the  summation  is  necessarily  non-negative,  the  strength  of  the  resulting 
spectrum  at  an  arbitrary  frequency  cannot  be  smaller  than  the  strength  of  the  weakest 
one  at  that  frequency. 

The  coherence  properties  of  the  source  are  manifested  by  the  distribution  of  the 
eigenvalues  (Xn  mfco)).  In  the  case  of  a  one-dimensional  Gaussian  Schell-model 
source,  for  example^. 


^  =  |pV2+l  +  p((P/2f  +  l]'^} (5.7) 


where  P  is  a  relative  correlation  length 


P  = 


0;(©)  ' 


(5.8) 
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Here  Oi((D)  is  the  r.m.s.  width  of  the  intensity  distribution  and  Oj^(a))  is  the  r.m.s. 
correlation  length  of  the  field  distribution  in  the  source  plane. 

In  this  chapter  we  only  consider  sources  which  contain  a  finite  number  of 
modes,  N.  Let  us  choose  a  source  with  eigenvalues  whose  frequency-dependence  is  of 
the  form 

=  (5.9) 

Here  s^^^co)  is  proportional  to  the  source  spectrum  S(x,  y,z-  0;©),  and  the  arc 
constants.  Explicitly, 

=  J S{x,r,<si)dxdy 

We  also  choose  the  modes  to  be  independent  of  the  frequency  ©,  i.e., 

=  %«.!«(»')•  (5.11) 

Since  the  domain  D  lies  in  a  plane  it  is  convenient  to  consider  modes  that  are  separable 
in  the  two  Cartesian  directions,  say  x  and  y,  i.c.  modes  of  the  form 

=  (5.12) 

For  the  functional  form  of  the  modes  we  choose  for  n  =  0,  1,  2... 
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where  //„  is  the  Hermite  polynomial  of  order  n  and  d,  is  a  constant  that  determines  the 
spatial  distribution  of  the  mode  in  the  j:-direction  [cf.  Ref.  5,  Eq.  (2.13)].  We  use  a 
stricdy  analogous  expression  for  <|)m(y). 

The  choice  of  frequency-independent  modes  isolates  the  frequency  dependent 
factors  in  the  expression  for  the  cross-spectral  density  of  the  field  in  the  source  plane; 
consequently,  for  an  arbitrary  source  spectrum  we  may  analyze  the  effect  of  the 
distribution  of  the  constants  ( A„)  on  the  resulting  field  spectrum.  This  is  equivalent  to 
examiiung  the  effects  of  spatial  correlations  of  the  source  fluctuations  on  the  spectrum 
of  the  field  that  the  source  produces. 

If  the  source  spectrum  and  the  field  spectrum  each  consist  of  a  single  line 
centered  at  frequencies  (Oq  and  0)’q  respectively,  one  may  characterize  the  spectral 
change  occurring  on  propagation  by  the  parameter 


Z  = 


0^0  ^*0 


CBr 


(5.14) 


However,  even  when  the  source  spectrum  consists  of  a  simple  line  such  as  a  Lorentzian 
or  a  Gaussian,  the  resulting  field  spectrum  is,  in  general,  a  distoned  line.  A 
comparison  between  the  spectral  lines  may  then  be  made  on  the  basis  of  peak 
frequencies  or  the  centroids  of  the  line.  In  the  following  sections  we  choose  the  later 
approach  for  calculating  the  new  center  fiwjucncy  ©’,  i.c.,  we  define  ©’  by  the  formula 


JcaS^^^tDldo) 
”  |5^“^(©)dQ) 


(5.15) 
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5.3  Far-zone  spectra 

The  far-zone  field  \tf„  „  of  a  typical  souru.  mode  Xo^n  can  be  obtained  by  substituting 
Eq.  (5.12)  into  Eq.  (5.6)  and  using  the  asymptotic  form  of  the  Green's  function.  We 
then  obtain 


,  .  InkcosQr  .r  ^ 

= - - - ^^(ku^y^JkUy), 


(5.17) 


where  6  is  the  angle  between  the  direction  of  observation  and  the  normal  to  the  source 
plane  and  ^„(/)  denotes  the  one-dimensional  spatial  Fourier  transform  of  the  function 
4>n(x), 

=  (5.18) 

On  substituting  Eq.  (5.13)  into  Eq.(5.18)  and  using  the  identity^ 


J  //,(a:)c"**^e^  dx  =  r^I2n  (5.19) 

we  obtain  the  following  expression  for  the  far  field  produced  by  the  mode  Xnjn- 
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(5.20) 


The  far-zone  spectrum‘s  of  the  field  is  now  obtained  by  substimting  Eq.  (5.9) 


into  Eq.  (5.5), 


(5.21) 


with  Yn  m  given  by  Eq.  (5.20).  We  may  now  express  the  formula  (5.21)  in  the 
following  form,  which  clearly  indicates  the  various  contributions  to  the  far-zone 
spectrum: 


2 

5(“\«;©)  =  5<®>(©)^^A/(«;d,;d^;©). 


(5.22) 


Equation  (5.22)  shows  that  the  far-zone  spectrum  is  a  product  of  the  source  spectrum,  a 
geometrical  factor,  and  a  “spectral  modifier”  A/(if ;  d,,  d^;  ©).  This  factor  is  given  by 


\  k  k  (  ku  \ 

M(u;d',d -,()})  =  — p=  ) - exp<-  -7=^=  +  -75^= 


XA  I  to  I 

"1^ 


(5.23) 
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To  examine  the  effect  of  the  spectral  modiner  on  the  spectrum  in  the  far  zone, 
let  us  consider  sources  whose  spectrum  consist  of  a  single  spectral  line,  centered  at 
frequency  (Dq.  We  define  the  relative  frequency  a  as  a  =  cc/coo,  and  the  characteristic 
mode  lengths  ^  ii  =  .  With  these  definitions  Eq.  (5.23)  may  be 

expressed  in  the  form 

exp|-a^j^(^,f  +  j| 

A  2 


For  a  fixed  direction  of  observation  the  spectral  modifier  M  is  seen  to  contain  an 
envelope  factor 


F  =  ^^exp|-a^  +(Tl«y)^  '  (5.25) 


and  a  weighted  sum  of  contributions  from  the  individual  modes.  As  we  indicated 
earlier  this  sum  does  not  involve  any  cross-mode  terms.  It  is,  therefore,  possible  to 
determine  the  spectral  changes  directly  from  the  knowledge  of  the  spectral  modifiers 
denoted  by  M,yn,  of  the  individual  modes; 


turn 


Ji2"'^'"n!m! 


exp|-a2|^(^u,f  + 


(5.26) 
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To  illustrate  the  behavior  of  the  spectral  modifier  of  individual  modes,  and  its 
effects  on  the  observed  spectra,  we  consider  sources  for  which  there  is  only  one-mode 
contribution  in  the  y-direction,  i.e., 

=  (5.27) 

We  also  limit  our  observation  points  to  the  plane  y  =  0,  where  the  spectral  modifier 
takes  the  form 

exp[-(a^J^].  (5.28) 

In  Figs.  5.1  and  5.2  we  show  examples  of  spectral  modifiers  of  single  modes 
for  selected  values  of  the  index  n  and  a  fixed  direction  of  observation.  One  can  see  by 
inspection  that  the  spectral  modifiers  of  modes  specified  by  n  =  4,  6,  9,  11,  give  rise 
to  blue  shifts  while  the  spectral  modifiers  for  modes  specified  by  n  =  5,  10,  12,  14 
give  rise  to  red  shifts. 
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Figure  5.1:  The  spectral  modifier  Af „  g  for  the  modes  n  =  4,  6,  9,  11  as  a  function  of  the 
relative  frequency  a.  The  positive  slope  of  the  curves  imply  a  blueshift  on  propagation. 


In  order  to  determine  the  type  of  spectral  shift  due  to  an  arbitrary  mode  n,,  we 
note  from  Eq.  (5.28)  that  the  fine  detail  in  the  spectral  modifier  is  due  to  the  square  of 
the  Hermite  polynomial.  The  behavior  of  one  such  factor  is  shown  in  Fig.  5.3.  We 


Figure  5.2:  Spectral  modifier  M,j^o  f'Q''  ^  modes  n  s  S,  10, 12, 14  as  a  function  of  the  relative 
frequency  a.  The  negative  slopes  of  the  curves  imply  a  ledshift  on  propagation. 
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Figure  5.3:  An  example  of  the  detail  in  the  spectral  modifier  due  to  the  square  of  the  Hermite 
polynomial  shown  here  for  n  =  1 1.  The  points  2^,,  denote  boundaries  of  regions  giving  rise 
to  particular  spectral  changes. 

note  that  at  the  center  frequency  a  =  1,  the  behavior  of  [//n(C)]^  reflects  its 
dependence  on  the  product  When  between  the  points  marked  by  the 

symbols  and  Cb»  diis  factor  produces  a  red  shift,  whereas  for  values  of  between 
and  ^  it  produces  a  blue  shift  When  the  values  of  correspond  to  the  point 
the  resulting  spectral  line  is  narrowed  while  for  values  of  ^  corresponding  to  the  point 
it  is  broadened.  It  is  clear  from  the  example  shown  in  Figs.  5.3  that  the  type  of 
spectral  effect  observed  for  a  fixed  value  of  is  determined  by  the  value  of 
compared  with  the  nearest  zero  of  the  Hermite  polynomial.  Specifically,  let  <  ^2 
two  consecutive  zeros  of  the  Hermite  polynomial,  and  let  be  the  value  at  which 
[//„(Q]2  attains  its  maximum  value  in  the  interval  ,  ^2)1  evidently,  ^2  • 

With  these  definitions  we  have 


Blue  shift 

Cl  <  Cm*<  Cm 

Red  shift 

Cm  <  C«^x<  C2 

line  splitting 

_ _ 

Table  5.1:  The  dqiendence  of  spectral  changes  on 
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Figure  5.4;  Normalized  far-zone  q)ectra  for  the  mode  n  =  7,  with  resulting  in  a 

blue-shifted  line  with  Z  =  -0.01 1. 


[Examples  of  the  three  cases  are  shown  in  Figs.  5.4-5.6.  The  source  qjectrum 
in  these  examples  was  taken  to  be  a  Lorentzian  line. 


Figure  5.5:  Normalized  far-zone  spectra  for  the  mode  n  =  7,  with  ^Xx  =  0.7,  resulting  in  a 
red-shifted  line  with  Z  =  0.0093. 
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Figure  5.6:  Normalized  far-zone  spectra  for  the  mode  n  =  7,  with  =  0.816  corresponding 
to  a  zero  of  the  Hermite  polynomial.  Line  splitting  is  evident  as  a  result  of  taking  close  to  the 
zero  of  H-j. 


/°^(a;5)  =  l/[l+(a-l)V5^].  (5.29) 


with  5  =  0.06. 


5.4  The  dependence  of  field  spectra  on  the  relative 
mode  strengths 

The  spectral  changes  which  occur  for  a  single  coherent  mode  is  a  consequence  of  wave 
propagation.  Such  effects  must  be  distinguished  from  the  spectral  changes  due  to 
source  correlations.  When  the  source  is  spatially  fully  coherent  it  consists  of  a  single 
coherent  tiKxle.  The  spectral  effects  which  then  arise  on  propagation  were  described  in 
the  previous  section.  However,  when  the  source  is  partially  coherent,  several  modes 
are  present  and  the  resulting  field  spectrum  then  also  depends  on  the  relative  strength  of 
the  modes  or,  equivalently,  on  the  distribution  of  the  constants  (A^). 
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Figure  5.7:  Spectral  modifiers  for  sources  characterized  by  the  spectral  modifiers 
M  =  ,  +  A/,  „  +  Af,  ,  +  Af„  „  (a).  aiKl  Af  =  Af^  ^  +  Af,„  „  +  Af,^„  +  Af,,  „  (b). 


Let  us  first  consider  a  source  consisting  of  two  coherent  modes  (ni.mj)  and 
(n2,m2).  It  follows  from  Eqs.  (5.22)  and  (5.26)  that  the  spectral  modifier  for  this 
source  is  just  the  sum  of  the  spectral  modifiers  of  the  two  modes,  i.e., 

A/(«;^,Ti;a)  =  ^(u;4,Ti;a)-«-A/^  ^(«;^,fi;a).  (5.30) 


In  terms  of  the  cross-spectral  density  Eq.  (5.1),  the  degree  of  correlation  of  the 
source  fluctuations  is  given  by 


p(r,,r2;(0) 


Vy^°>(r,,r2;Q)) 

-^VV(°>(r,.r,;(D)^H^^®^(r2.r2;©)' 


(5.31) 
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In  Fig.  5.7  we  show  examples  of  spectral  modifiers  for  sources  consisting  of 
four  modes.  The  degree  of  correlation  p.f'-j,  r2;  to)  for  the  <  orresponding  sources  is 
shown  in  Figs.  5.8  and  5.9  for  points  along  the  x-axis,  with  Tj  =  (C,  0)  and 


0.0  1.0  2.0  3.0  4.0  5.0 

Spatial  offset 


Fig.  5.8;  Degree  of  spatial  correlation  p(0,  5i)  of  the  two  sources  shown  in  Fig.  5.7  (a), 
where  the  spatial  offset  is  &  =  x^j2d^  . 


0.0  1.0  2.0  3.0  4.0  5.0 

Spatial  offset 


Fig.  5.9;  Degree  of  spatial  correlation  u(0.  &)  of  the  two  sources  shown  in  Fig.  5.7  (b), 
where  the  spatial  offset  is  fix  =  . 
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Figure  5.10:  Comparison  of  the  ^lectral  modifiers  for  a  source  consisting  of  modes  0-30  (a) 
and  a  source  consisting  of  the  single  mode  [n  =  9,  m  =  0]  (b).  The  normalization  ensures  that 
both  modifiers  have  the  same  magnitude  at  the  center  frequency. 


r2  =  (x,  0).  In  some  cases  the  degree  of  correlation  may  approach  a  constant  value 
with  increasing  separation  between  the  two  source  points.  This  result  is  well  known 
for  Gaussian  Schell-model  beams*. 

When  the  source  consists  of  many  modes  taken  with  equal  weights,  it  is 
equivalent  to  a  spatially  incoherent  source  [cf.  ref.  5].  Figure  5.10  shows  the  spectral 
modifier  for  a  source  consisting  of  modes  n  =  0  through  n  =  30.  For  comparison  we 
also  show  the  spectral  modifier  for  a  source  consisting  of  the  single  mode  n  =  9.  It 
can  be  seen  from  this  figure  that  when  the  source  consists  of  many  modes,  the  spectral 
modifier  is  relatively  constant  and  consequently  the  far-zone  spectrum  produced  by  this 
source  cannot  differ  appreciably  finom  the  source  spectrum. 


5.5  Spectral  changes  in  the  near  zone 

The  spectrum  of  the  field  throughout  the  half-space  z  >  0  may  be  evaluated  by  using 
the  angular  spectrum  representation  for  each  one  of  the  modes,  i.e.. 
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V„.;n 


(r-.co)  =  JK  ^(/>.<7;£D)c 


ik{px+qy+mz) 


dpdq. 


(5.32) 


In  Eq.  (5.32)  r  =  ix,y,z>  0), 

j  J (5.33) 


and 


when 

(a) 

when 

+q^  ^1 . 

ib) 

(5.34) 


It  can  be  shown  by  numerical  evaluation  that  when  ib  >  100  one  may  neglect, 
to  a  good  approximation,  the  contributions  from  evanescent  waves,  i.e.  waves  for 
which  +  ^2  >  1  Tjie  infinite  domain  of  integration  in  Eq.  (5.32)  can  then  be 
replaced  by  the  finite  domain  iP-  +  q^<\. 

To  simplifj'  the  calculation  we  consider  the  source  modes  given  by  Eq.(5.12)  in 
Eq.  (5.33).  The  integration  over  the  domain  D  may  be  carried  out  in  closed  form  by 
extending  the  domain  of  integration  over  the  whole  x,y*plane  and  using  the  identity 
(5.19).  One  then  finds  that 
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\JPr<I>CO)  =  ~ 


(5.35) 


The  fields  \yn,in(r;co)  are  obtained  by  substituting  from  Eq.  (5.35)  in  Eq.  (5.32)  and 
letting  dj^  =  dy  =  d.  The  result  is: 


V^«(x.y.z;o)) 


2jiJ  V  d  -^2'‘n! 


(5.36) 


p^+q^il 


The  integral  (5.36)  has  been  evaluated  by  many  authors  within  the  accuracy  of 
the  paraxial  approximation^.  It  has  also  been  evaluated  under  noore  general  conditions 
or  by  making  use  of  cylindrical  symmetry^®  Because  we  are  interested  in  observation 
points  that  may  lie  outside  the  paraxial  region,  we  evaluated  the  integrals  (5.5) 
numerically  for  selected  values  of  the  mode  index  n. 

In  Table  5.2  we  compare  the  relative  frequency  shifts  calculated  for  the  near 
zone  and  for  the  far  zone. 
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Mode 

Near  zone 
Z-shift 

Far  zone 
center 
fieouenev 

Far  zone 
Z-shift 

5,0 

0.992737 

-0.007316 

0.992654 

-0.007399 

0.000083 

9,0 

1.020491 

0.020079 

1.021117 

0.020681 

0.000602 

11,0 

1.011317 

0.01119 

1.011482 

0.011352 

0.000162 

12,0 

0.993354 

-0.00669 

0.993168 

-0.006879 

0.000189 

13,0 

1.007263 

0.007211 

1.007341 

0.007288 

0.000077 

14,0 

0.985337 

-0.014881 

0-9845217 

-0.015722 

0.000841 

15,0 

1.004433 

0.004414 

1.0044527 

0.004433 

0.000019 

Table  5.2:  The  Z-numbers  (relative  frequency  shifts)  for  observation  points  in  the  far  and  in 

the  near  zone,  for  selected  modes. 


The  rightmost  column  in  Table  5.2  indicates  that,  at  least  for  the  modes  considered 
here,  most  of  the  spectral  shift  is  already  present  for  observation  points  in  the  near 
zone.  This  observation  is  confirmed  in  Chapter  6. 


5.6  Summary 

In  this  ch^ter  we  employed  the  coherent-mode  representation  to  analyze  the  changes  in 
the  spectrum  of  light  which  is  generated  by  a  class  of  secondary,  partially  coherent, 
planar,  secondary  sources.  Our  treatment  isolated  the  spectral  effects  due  to  the 
individual  modes  and  has  demonstrated  that  even  coherent  fields  exhibit  a  non- 
negligible  spectral  shifts  for  some  directions  of  observation.  These  may  have  to  be 
taken  into  account  in  various  applications,  for  example,  in  determining  the  speed  of 
moving  objects  from  reflected  light  on  the  basis  of  spectral  line  shifts. 
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The  effect  of  the  state  of  coherence  of  the  source  on  the  observed  field  spectra 
was  examined  by  varying  the  eigenvalues  A„  We  found  that  when  the  source 
consists  of  many  coherent  modes  with  the  same  eigenvalues  (corresponding  to  the 
incoherent  limit),  the  spectrum  of  the  field  does  not  differ  appreciably  from  that  of  the 
source. 

The  changes  in  the  spectrum  of  the  field  on  propagation  from  the  immediate 
neighborhood  of  the  source  plane  to  the  far  zone  were  studied  by  calculations  of  the 
spectrum  for  distance  kz  =  100  and  for  kz  — >  «>o.  We  found  that  in  the  case  we 
studied,  the  spectrum  close  to  the  source  {kz  =  100)  already  possesses  most  of  the 
features  that  are  present  in  the  spectrum  of  the  field  in  the  far  zone. 

The  evaluation  of  the  near  zone  spectrum  in  Sec.  5.5  is  computationally 
intensive.  In  Chapter  6  we  use  the  paraxial  approximation  to  perform  this  calculation. 
The  paraxial  approximation  significantly  reduces  the  difficulty  involved  in  the 
computation  which  gives  a  genuine  insight  into  the  development  of  spectral  changes 
with  propagation  distance. 
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6.1  Introduction 

The  work  described  in  the  previous  chapters  has  focused  on  light  propagation  in  free 
space  O’  in  a  rarefied  scattering  medium.  In  this  chapter  we  examine  the  changes  in  the 
spectrum  of  the  field  occurring  on  propagation  through  homogeneous  and 
inhomogeneous  dispersive  media.  In  particular,  we  investigate  the  development  of  the 
spectrum  as  a  function  of  the  propagation  distance  and  its  dependence  on  the  state  of 
coherence  of  the  source  in  both  homogeneous  and  inhomogeneous  media.  The 
inhomogeneous  medium  considered  here  is  a  graded-index  medium  whose  refractive 
index  varies  quadratically  in  the  radial  direction^  Such  a  medium  is  readily  available 
for  experiments  in  the  form  of  so-called  Selfoc  fibers. 

The  changes  in  the  coherence  properties  of  light  propagating  through  various 
types  of  waveguides  have  been  investigated  by  many  authors^.  In  most  of  the  work 
encountered  in  the  literature,  the  state  of  coherence  is  characterized  by  the  mutual 
coherence  function.  As  we  explained  in  Chapter  1,  this  approach  is  not  very  suitable 
for  the  examination  of  spectral  changes.  Agrawal  et  al.  [  Ref.  2  (a)]  considered  how 
the  cross-spectral  density  of  the  incident  light  changes  on  propagation  in  such 
multimode  fibers.  We  use  the  analysis  of  that  reference  to  derive  a  closed-form 
expression  for  the  spectrum  of  the  field  at  an  arlntrary  distance  from  the  source. 

The  general  expression  derived  for  a  graded-index  medium  can  be  used  to 
analyze  the  spectral  changes  occurring  in  a  dispersive  homogeneous  medium  in  the 
appropriate  limit  The  later  result  reduces  to  the  well  known  free-space  result^  in  the 
limit  in  which  the  refractive  index  is  unity.  Our  expression  is,  however,  valid  for 
arbitrary  propagation  distances  and  allows  us  to  examine  how  the  spectrum  evolves 
from  the  near-field  to  the  far-field  region.  We  illustrate  our  results  by  using  physical 
parameters  that  apply  in  many  practical  configurations.  In  particular,  we  show  how  the 
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spectrum  of  light  can  be  shifted  toward  the  shorter  or  the  longer  wavelength,  depending 
on  the  propagation  distance  and  the  state  of  coherence  of  the  source.  Our  results 
indicate  that  the  spectral  shift  occurring  in  the  far-zone  region  is  considerably  enhanced 
in  a  homogeneous  medium,  with  the  enhancement  factor  depending  on  the  index  of 
refraction  of  the  medium. 

We  examine  the  evolution  of  spectral  shifts  from  the  near  to  the  far  zone  by 
deriving  an  alternative  expression  for  the  field  spectrum  that  is  valid,  within  the  paraxial 
approximation,  for  an  arbitrary  propagation  distance  in  free  space.  We  use  this 
expression  to  obtain  the  spectral  shift  for  optical  fields  generated  by  a  Gaussian  Schell- 
model  source,  and  to  study  how  the  shift  changes  during  transition  from  the  near  to  the 
far  zone. 


6.2  Propagation  of  the  spectrum  in  graded  index  fibers 

In  this  section  we  derive  expressions  for  the  propagation  of  the  cross-spectral  density 
of  the  field  in  graded-index  fibers.  We  establish  the  notation  and  derive  an  expression 
for  the  field  spectrum. 

Qxisider  a  graded-index  fiber  with  the  axis  of  symmetry  along  the  z-direction 
(see  Fig.  6.1).  The  fiber  is  characterized  by  an  index  of  refraction  having  the  parabolic 
profile 


n^(x,r,(0)  = 


/iJ(to)[l-a^(tu)(x^-»-y^)] 


for  +  (a) 

for  x^+y^>/^.  (b) 


(6.1) 
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4 

Figure  6.1:  Illustrating  the  geometry  and  the  notation.A  point  in  the  source  plane  z  =  0  is 
denoted  by  (^,  t))  and  an  obsovadon  point  is  denoted  by  (x,  y,  z). 


where  Rq  is  the  core  radius,  Hq  is  the  index  of  refraction  at  the  center  of  the  fiber,  a  is 
the  radial  gradient  of  the  index,  and  (0  =  kc  (c  =  vacuum  speed  of  light)  is  the 
frequency  associated  with  the  free-space  wavenumber  k. 

To  obtain  an  expression  for  field  propagation  in  this  medium^  we  first  make  the 
assumption  that  Eq.  (6.1a)  is  valid  for  all  x  and  y.  It  follows  that  the  differential 
equation  governing  the  propagation  of  a  field  V  in  this  fiber  is  in  the  form 

{v2+ifeJ[l-a2(jc2+/)]jv(x,y,z)  =  0  .  (6.2) 
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Its  solution  can  be  expressed  as  an  infmite  series  of  Hermite  Gaussian  functions^ 


V(x,y,2)=  X 

n,in=0 


rvjy'l 

Uof" 

1^0  J 

’  (6.3) 


where  Wq  is  the  spot  size 


WQ  =  ^2/ka, 

and  are  the  propagation  constants  given  by 

Pn.m  =  ^k^-2ka{n  +  m  +  \). 


(6.4) 


(6.5) 


In  Eq.  (6.3)  the  coefficients  are  determined  from  the  boundary  value  of  the  field 
incident  on  the  fiber  0)  using  the  orthogonality  of  the  Hennite-Gaussian 

functicais.  One  then  finds  that 


n,m 


2  1 


xJv(^.Ti,0)«, 


H. 


yv, 


0  y 


exp 


Wn 


(6.6) 


d^dr]. 


where  ^  and  ti  are  the  coordinates  in  the  input  plane  D. 

If  follows  from  Eqs.  (6.3)  and  (6.6)  that  the  field  V(r)  after  propagation  to  any 
plane  z  >  0  may  be  expressed  in  the  form 
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V(r)  =  j  V{^,n,0)K{r-^,r])d^  dr],  (6.7) 

D 

where  the  propagation  kernel  K(r,^,T])  is  given  by 

K(r;^,Tj)  =  ^e‘^^'^^^cxpl-^^l^cosazU^  +  T]^]-{^  +  T]y)  |  ,  (6.8) 

2ni  sin  az  (sinazL2  ^  ’  JJ 

and 

(p(r)  =  k\z  +  — — - — cosazfx^+y^]  .  (6.9) 

L  2sinaz  ' 

When  the  medium  is  homogeneous,  a  =  0  and  Eq.  (6.8)  reduces  to 

which  is  also  the  propagation  kernel  for  the  paraxial  Fresnel  approximation. 

Using  Eqs.  (6.7)  and  (6.8),  the  cross-spectral  density  of  the  field  for  any  two 
points  in  the  fiber  is  given  by 

W(ri.r2;co)  =  JJ  A'*(ri,fl;0))/s:(r2,ft;©)W(/^./^;CD)d^^  d'^p^ ,  (6.1 1) 

D 
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where  Vy(p^,p^;to)  is  the  cross-spectral  density  in  the  source  plane  and  Pj  =  is 

a  radius  vector  in  that  plane.  The  spectrum^  of  the  field  is  obtained  by  setting 
ri  =  r2  =  r  in  Eq.  (6.11),  i.e., 

S(r;a))  =  JJ K\r,p-(o)Kir,p^;(ii)W(p^,p^;ai)d'^p^  d'^p^.  (6 j2) 


On  substituting  from  Eqs.  (6.8)  and  (6.9)  into  Eq.  (6.12)  we  obtain  the  following 
expression  for  the  spectrum  of  the  field: 


ika 

'cosotz/t; 

Isinotz 

.  2 

(6.13) 


Equation  (6.13)  can  be  used  to  obtain  the  field  spectrum  at  any  point  r  for  a 
given  form  of  the  cross-spectral  density  of  the  incident  field  in  the  source  i.  ’me  z  =  0. 
We  will  choose  a  Gaussian  Schell-model  for  the  cross- spectral  density  in  the  source 
plane,  viz.. 


W(Pj,p!^;a))  =  S^°^(o))ftxp 


xcxp 


(6.14) 
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where  is  the  r.m.s  width  of  the  Gaussian  intensity  distribution  [Full-Width-Half- 
Maximum(FWHM)  =  2.35  Oj]  and  <Tg  is  the  r.m.s  width  of  the  spatial  correlation. 
On  substituting  Eq.  (6.13)  into  Eq.  (6.14)  and  performing  the  integration  [see 
Appendix  (T]  we  obtain  the  following  expression  for  the  spectrum  of  the  field: 

5(r;o))  =  S^®^(o))M(r;a;ti)),  (6.15) 


where  the  si>ectral  modifier  M  is  given  by 


A/(r;a;co): 


(  to,  f 

1  A  J  ’’ 

(6.16) 


In  Eq.  (6.16)  we  have  used  the  notation 

A  =  2abO;Sin(xz/a2  , 


(6.17) 

(6.18) 


and 


tocoscczV 
2asinaz  ) 


(6.19) 


Equations  (6.15)  and  (6.16)  are  general  expressions,  valid  within  the  paraxial 
approximation,  for  the  spectrum  of  the  field  produced  by  planar  secondary  Gaussian 
Schell-model  sources.  In  particular,  Eq.  (6.15)  is  valid  at  any  distance  z  from  the 
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source  plane  and  for  any  range  of  real  valued  n(a))  and  ce(co).  It  can  be  easily  verified 
that  5(r;  to)  reduces  to  S^^H(o)cxp^-r^/2ajj  in  the  limit  z  ->  0,  as  one  might 
expect.  In  the  limit  when  a  0  our  formulation  corresponds  to  the  case  of 
dispersive  homogene  )us  media.  The  spectral  modifier  for  homogeneous  media  My,  is 
given  by 


My,(r,co) 


yexd 


+  1 


(6.20) 


where 

z  2/fca^ 


(6.21) 


The  parameter  Zj  is  analogous  to  the  so-called  diffraction  length  or  the  Rayleigh  range^ 
encountered  in  propagation  of  coherent  Gaussian  beams  and  reduces  to  it  in  the 
coherent  limit  (<jg  »  q). 

Equation  (6.20)  is  also  applicable  to  paraxial  free-space  propagation  if  we  set 
/i((o)  =  1;  then  k-k^.  Since  we  have  examined  in  the  previous  chapter  the  spectral 
changes  that  occur  on  free-space  propagation,  we  consider  this  case  first  Later  we  will 
consider  the  general  case  in  which  nfto)  >  1  and  a  >  0.  In  order  to  make  our  results 
readily  available  for  exp)erimental  verification,  we  assume  that  the  source  spectrum 
S(0)(©)  corresponds  to  that  of  a  Gallium  Phosphide  (GaP)  visible  light  source*,  which 
is  well  approximated  by  a  Lorentzian  line  centered  at  564  nm 
(Vq  =  C0q/2ti  =  532  THz)  and  having  a  FWHM  of  36  nm  (=  34  THz).  All  the 
numerical  results  in  this  chapter  pertain  to  this  source  spectrum. 
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6.3  Free-space  propagation 

In  the  previous  chapters  we  examined  the  changes  in  the  spectrum  of  light  on 
propagation  in  free  space  and  we  reviewed  the  relevant  published  work.  So  far  there 
has  not  been  a  clear  understanding  of  the  development  of  the  spectral  changes  with 
increasing  propagation  distance  from  the  source,  a  fact  that  may  be  due  to  the 
complexity  of  the  computations  involved^.  Our  analysis  provides  a  relatively  simple 
way  for  understanding  the  transition  from  near  to  far  field  within  the  paraxial 
approximation. 

In  our  notation  the  spectrum  of  the  field  after  propagating  a  distance  z  in  free- 
space  is  given  by 


5^(r;o))  =  S(°>(a))M^(r;tD)  ,  (6.22) 


where  is  the  spectral  modifier  for  free-space  propagation  that  is  obtained  from  of 
Eq,  (6.20)  by  setting  k  = 

The  expression  for  the  spectrum  of  the  field  in  the  far- zone  is  obtained  in  the 
limit  koz  — >  oo  (with  fixed  direction  of  observation),  and  it  is  given  by 


5j-~\r,©) 


exp 


4a" 


(6.23) 


where  the  superscript  («»)  indicates  the  far-zone  limit.  Equation  (6.23)  is  in  agreement, 
within  the  paraxial  approximation,  with  a  known  result  fw  far-zone  radiant  intensity  of 
Gaussian  Schell-model  sources^. 
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Figure  6.3:  Normalized  spectral  modifler  Aff  for  propagation  distance  Aqz  =  100  in  free 
space.  The  spectral  modifier  is  shown  as  a  function  of  frequency  v  for  =  20  and  four  different 
values  of  the  ctMrelation  length;  kQOg  =  1.0  (a).  kQOg  =  8.0  (b),  k^Og  =  10  (c)  and  =  20 
(d).  The  direction  of  the  spectral  shift  is  determined  by  the  slope  of  at  the  center  frequency  of 
the  source. 


We  now  return  to  the  general  expression  for  the  spectrum  of  the  field  in  free 
space  [Eq.  (6.20)],  and  evaluate  the  spectral  modifier  My  for  sources  with  different 
states  of  coherence,  governed  by  the  values  of  Og  and  Oj. 

Figure  6.3  shows  the  variation  of  the  spectral  modifier  My  with  the  frequency 
V  =  ©/27t  for  several  choices  of  <Tg  and  <7j,  when  the  propagation  distance  is 
z  =  1000A,o/2n  (i.e.,  when  koz  =  1000),  and  the  observation  direction  makes  an 
angle  of  10*  with  the  z-axis.  Each  one  of  the  curves  is  normalized  so  that  its  maximum 
value  is  unity.  It  follows  from  Eq.  (6.20)  that  when  the  specL*al  modifier,  considered 
as  a  function  of  v,  has  a  positive  slope  at  v  =  Vq,  the  resulting  line  is  blue-shifted 
while  a  negative  slope  of  the  spectral  modifier  at  that  frequency  results  is  a  red-shifted 
spectrum.  If  the  spectral  modifier  is  not  uniform  throughout  the  frequency  range  of  the 
source  spectrum,  the  nature  of  the  spectral  changes  may  be  more  complicated  (for  more 
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Figure  6.4:  Normalized  spectral  modifier  Mf  for  (Hopagadcn  distances  io:  -  100  (a),  tor  =  250 
(b),  and  jbQZ  =  600  (c)  in  free  space  for  tot^i  =  20,  =  10-  '’0  =  532  THz  a  blue  shift  is 

obtained  ftv  =  100  and  a  red  shift  for  kQZ  =  600. 

details  sec  Sec.  5.3).  The  examples  shown  in  Fig.  6.3  correspond  to  a  blue  shift  for 
ikotTg  <  9  and  red  shift  for  k  ^  ^  10  for  GaP  source  for  which  Vq  =  532  THz. 

For  a  given  state  of  coherence  the  spectral  shift  also  depends  on  the  propagation 
distance.  This  is  shown  in  Fig.  6.4  where  wc  compare  the  spectral  modifier  for 
different  propagation  distances  from  sources  with  ifcoOj  =  20  and  =10.  The 
spectral  shift  is  to  vaid  higher  frequencies  for  kQZ  =100  (positive  slope  at  n  =  /iq), 
and  toward  lower  frequencies  for  kQZ  =  600  (negative  slope  at  n  =  /Iq).  Figure  6.5 
shows  the  spectra  of  the  GaP  source  when  =  10  (i.c.  in  the  near  zone)  and  when 
koz  =  1000  (in  the  far  zone)  with  kQOi  =  20  and  kQO^  =  20.  The  source  spectrum 
exhibits  a  red  shift  in  the  far  zone,  and  a  slight  blue  shift  in  the  near  zone.  It  should 
also  be  noted  that  the  spectrum  becomes  asymmetric  as  a  result  of  propagation. 

We  quantify  the  magnitude  of  the  spectral  shift  by  defining  a  parameter  Av  that 
corresponds  to  the  shift  of  the  spectral  peak  from  from  the  location  of  the  peak  of  the 
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Fignre  6.S:  Nonnali^cd  field  spectrum  for  observation  at  an  angle  of  10’  off  axis  and  a 
propagation  distance  kgz  =  1000.  The  source  is  characterized  by  iloa]  s  20  and  ito^g  ~ 
solid  line  shows  the  original  source  spectrum  and  the  dashed  line  shows  the  red-shifted  field  spectrum. 

source  spectrum.  In  Fig.  6.6  we  show  the  spectral  shift  Av  as  a  function  of  the 
propagation  distance  koz,  for  koSj  =  20  and  several  values  of  kgOg.  We  note  that 
when  the  source  is  relatively  incoherent  (koOg  ^  1)  the  spectral  shift,  which  is 
towards  the  blue  for  this  angle  of  observation,  develops  rapidly  with  propagation 
distance.  When  the  source  is  relatively  coherent  (^<Tg  »  1)  an  initial  blue  shift  turns 
into  a  red  shift  with  increasing  koZ. 

For  the  states  of  coherence  considered  in  this  paper,  the  transition  from  the  blue 
shift  to  the  red  shift  takes  place  when  koZ  ~  100.  In  all  cases  the  frequency  shift  in  the 
far  zone  approaches  an  asymptotic  value  which  depends  on  The  constant  value  . 
The  frequency  shift  for  kgOg  =  25  [curve  (d)  in  Fig.  6.6]  is  about  10%  of  the  source 
spectral  width  (FWHM  =  34  THz). 
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Propagation  Distance  k^z 

Figure  6.6:  Frequency  shifts  Av  versus  propagation  distance  fcM'  sources  chr^terized  by  the 
same  value  of  kQOi  =  20  and  different  values  of  kQOgi  kQOa  =  1  (a),  =  10  (b),  )^)a»  =  20  (c; 

and  iko<Jg  =  25  (d). 


6.4.  Spectral  changes  in  homogeneous  media 
In  the  previous  section  we  showed  that  the  changes  in  the  spectrum  on  propagation  in 
free  space  depend  on  the  state  of  coherence  and  on  the  propagation  distance  from  the 
source.  In  this  section  we  consider  propagation  through  homogeneous  media  for 
which  the  index  of  refraction  [nfco)  >  1]  is  independent  of  position  in  space.  We 
assume  that  the  index  of  refraction  is  frequency  dependent,  a  feature  that  indicates  the 
dispersive  nature  of  the  homogeneous  medium.  The  wave  number  is  then  given  by 

k  =  «(©)-.  (6.24) 

c 
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On  substituting  Eq.  (6.24)  in  Eq.  (6.20)  we  obtain  the  expression  for  the 
spectral  modifier  for  dispersive  homogeneous  media.  When  the  refractive  index  n(a))  is 
nearly  constant  over  the  source  spectral  width,  the  medium  acts  as  a  non-dispersive 
homogeneous  medium  of  constant  refractive  index  hq  =  n(a\)),  where  coo  is  the  central 
frequency  of  the  source  spectrum.  We  will  consider  the  non-dispersive  case  first 

It  is  evident  from  Eq.  (6.20)  tliat  the  spectral  modifier  Mj,  for  a  non-dispersive 
homogeneous  medium  is  identical  to  that  of  free  space  if  /cq  is  replaced  by  ngko-  Thus, 
the  free  space  results  for  the  spectral  modifier  shown  in  Figs.  6.3  and  6.4  apply, 
provided  that  the  scaling  factor  kg  is  appropriately  modified.  The  spectral  changes  can 
be  quite  different  as  a  result  of  the  scaling.  The  comparison  between  the  spectral 
changes  occurring  on  propagation  in  free-space  and  on  propagation  in  homogeneous 
non-dispersive  media  must  be  considered  separately  for  distances  short  and  long 
compared  with  z^.  For  a  short  propagation  range  (z  «  z^) 


1 


1 . 


(6.25) 


and  the  spectral  modifier  Mf,  becomes  independent  of  hq.  Hence  for  such  propagation 
distances  we  expect  no  difference  between  the  spectral  changes  occurring  in  dispersive 
homogeneous  media  and  free  space.  For  long  propagation  distances  (z  »  Zj)  the 
approximation 


(6.26) 
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holds,  and  the  spectral  modifier  for  propagation  in  homogeneous  media  [n(CD)  >  1] 
then  differs  by  a  multiplicative  factor 


2  2 
x+y 


(6.27) 


in  comparison  with  the  spectral  modifier  for  fiee-space  propagation.  In  Eq.  (6.27)  zq  is 
given  by  the  expression 


(6.28) 


The  extent  of  the  spectral  changes  taking  place  under  these  circumstances  depends  on 
the  value  of  n©.  Figure  6.7  gives  a  comparison  of  spectral  shifts  in  free  space  (wq  =  1) 
and  in  two  non-dispersive  homogeneous  media  of  refractive  indices  =  1.5  and  2.0 
when  k^Ci  -  20  and  k^cf^  =10.  The  most  notable  featiue  is  that  the  far- zone  value 
of  the  spectral  shift  increases  with  increasing  refractive  index  hq.  This  is  an  important 
feature  which  shows  that  the  spectral  changes  are  enhanced  in  a  homogeneous  medium. 


Spectral  changes  on  propagation  of  partially  coherent  light  in  media 


Chap(er6 


157 


Figure  6.7:  Comparison  of  frequency  shifts  for  propagation  in  non-dispersive  homogeneous  media. 
The  frequency  shifts  for  a  fixed  angle  of  observation  (10*)  are  shown  for  prc^gation  in  free  space  (a),  for 
propagation  in  a  homogeneous  medium  of  an  index  of  refraction  nftuo)  =  1.5  (b)  and  for  propagation  in  a 
medium  of  index  of  refraction  /i((ao)  =  2.0  (c).  The  observation  angle  is  10*  and  the  source  parameters 
are  k^Oi  =  20  and  ibgOg  =  10. 


We  next  consider  propagation  in  dispersive  homogeneous  media.  As  in  the 
case  of  non-dispersive  media,  we  must  consider  the  changes  in  the  spectrum  for  short 
and  long  propagation  distances  separately.  For  short  propagation  distances  (z  «  z^), 
the  quantity  z^j  z^  jis  again  independent  of  the  index  of  refraction  and  the  spectral 

effects  arc  identical  to  those  encountered  in  frec-space  propagation.  For  long 
propagation  distances  (z  »  zq). 


/(©)  =  nJ(CD)cxpj 


(6.29) 


and  hence  the  difference  in  the  spectral  effects  in  this  medium  from  those  generated  in 
free-space  depend  on  the  variation  of  nfeo)  in  the  frequency  range  covered  by  the  source 
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Figure  6.8:  Comparison  of  frequency  shifts  for  di^rsive  homogeneous  media,  is  shown  as 
a  finction  of  kQZ  for  propagation  in  free  space  (a),  propagation  in  pure  silica  (b),  and  propagation  in 
silica  doped  with  7.9%  Ge02  (c).  The  observation  angle  is  10*  and  the  source  parameters  are 
igOj  =  20  and  kQffg  =  10  . 


spectrum  [cf.  Eq.  (6.27)].  We  illustrate  our  results  by  using  a  slab  of  silica  glass  as  an 
example  of  dispersive  homogeneous  medium. 

Figure  6.8  shows  the  frequency  shift  An  obtained  after  light  from  GaP  source 
propagates  through  a  slab  of  silica  glass  of  various  thicknesses.  The  source  parameters 
are  kgOi  =  20  and  =  10.  The  frequency  dependence  of  n((o)  was  obtained  by 
using  the  well-known  Sellmeier  formulaJO 


n  (0))-l  +  V— y 


B(o^ 

J  J 


“  to  -  to 

j=\  ‘"y 


(6.30) 


For  pure  silica  B j  =  0.6961663,  B j  =  0-^079426,  B3  =  0.8974794, 
Aj  =  0.0684043,  A2  =  0.1 162414,  A3  =  9.896161  mm,  where  Aj  =27tc/a)j. 

The  effect  of  dopants  on  the  spectral  shift  can  be  easily  included  in  our  analysis. 
For  example,  the  refractive  index  n(to)  of  silica  glass  can  be  increased  by  doping  it  with 
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germania  (GeC)2)-  The  refractive  index  n(co)  is  still  given  by  the  Sellmeier  formula  but 
the  parameters  and  cOj  are  different  and  depend  on  the  amount  of  the  dopant  As  an 
exanqjle,  we  consider  silica  glass  doped  with  7.9%  Ge02,  fOT  which  the  parameters  are 
=  0.7136824,  B2  =  0-4254807,  B3  =  0.8964226,  Aj  =  0.0617 167, 
A2  =  0.1270814,  A3  =9.896161  mm.  Figure  6.8  shows  the  expected  change 
(dashed  curve)  in  the  fiequency  shift  The  shift  is  slightly  larger  for  doped  silica  since 
the  dopant  increases  the  refractive  index  by  a  small  amount  In  both  cases  (pure  silica 
and  slightly  doped  silica)  the  frequency  shifts  in  tfte  far  zone  are  much  larger  than  those 
that  would  be  produced  in  free  space. 

The  main  conclusion  of  this  section  is  that  correlation-induced  spectral  shifts  are 
enhanced  in  a  homogeneous  medium  of  refractive  index  n(Q))  >  1.  The  fiequency 
dependence  of  the  reft^tive  index  n(o))  is  not  critical  as  the  enhancement  is  found  to 
occur  even  when  n  is  frequency  independent  The  origin  of  the  enhancement  factor  for 
large  propagation  distances  can  be  understood  by  referring  to  Eq.  (6.20)  and  using  the 
formula  k  =  a)n(0))/c  in  Eq.  (6.21).  The  Gaussian  factor  in  Eq.  (6.20),  plotted  as  a 
function  of  co,  is  narrower  for  a  homogeneous  medium  than  for  free  space.  It  is  this 
feature  of  the  spectral  modifier  that  is  responsible  for  a  larger  spectral  shift  when 

/l(CD)  >  1. 


6.5  Propagation  in  inhomogeneous  media 

In  this  section  we  return  to  the  general  expression  for  the  spectrum  of  light  in  a  graded- 
index  fiber  [Eqs.  (6.15)-(6.19)].  In  this  case  the  parameter  a  is  non- zero;  it  depends 
on  the  fiber  design.  In  particular  we  consider  a  fiber  whose  core  is  made  of  doped 
silica  (7.9%  Ge02  at  the  core  center)  and  a  cladding  made  of  pure  Si02.  If  /ii(o))  is  the 
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Figure  6.9:  Frequency  shift  Av  versus  the  propagation  distaiKe  Ji()Z  in  a  dispersive  graded-index 
medium  (a).  Curve  (b)  shows  Av  when  the  inhomogeneous  nature  of  the  medium  is  ignored  by  settir.g 
a  =  0.  Curve  (c)  shows  the  corresponding  result  for  Cree-space  propagation.  The  observation  angle  is 
10*  and  the  source  parameters  were  chosen  to  be  =  20  and  =  10  . 


refractive  index  at  the  core  center  (r  =  0)  and  n2((o)  is  the  refractive  index  at  the 
boundary  (r  =  Rq),  the  parameter  a  is  given  by 

a(co)  =  -^-^l-n|(a))/nf(to)  .  (6.31) 


Since  /ti(a))  and  /i2(CD)  can  be  obtained  by  using  Eq.  (6.29),  the  frequency  dependence 
of  oc(ci))  is  readily  determined.  In  the  following  calculations  we  take  the  core  radius  to 
be  i?o  =  25  pm. 

Figure  6.9  shows  the  frequency  shift  as  a  function  of  the  propagation  distance 
for  a  0  in  such  a  medium,  and  the  shifts  generated  in  a  homogeneous  medium  of 
refractive  index  n].  We  note  that  in  the  range  of  propagation  distances  shown  in  the 
figure,  the  frequency  shift  is  larger  for  a  graded-index  medium  than  in  the 
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homogeneous  medium.  This  enhan''ement  of  the  frequency  shift  is  due  to  the 
inhomogeneous  nature  of  the  medium  and  depends  on  the  functional  form  of  the 
inhomogeneity. 

Although  the  frequency  shift  A\  tends  to  a  constant  (far-zone)  in  a 
homogeneous  medium,  no  such  limit  exists  in  the  graded-index  medium  considered 
here.  This  feature  is  due  to  the  imaging  property  of  a  medium  with  a  quadratic  variation 
of  the  refractive  index.  Such  a  medium  reproduces  the  incident  field  periodically,  with 
a  period  given  by  Az  -  2n/a,  a  feature  that  is  due  to  the  periodic  nature  of  the 
propagation  kernel,  Eq.  (6.8).  One  would  thus  expect  that  the  spectrum  S(r;  co)  given 
by  Eq.  (6.12)  also  reduces  to  the  source  spectrum  for  z  =  Irrnila,  where  m  is  a 
positive  integer.  We  show  in  Appendix  D  that  this  is  indeed  the  case.  Furthermore,  we 
fmd  that  the  source  spectrum  is  reproduced  not  only  at  z  =  Imitla,  but  also  at 
z  =  (2m  +  1  )n/a,  except  for  a  spatial  inversion  of  the  intensity  distribution.  For  a 
symmetric  intensity  profile  such  as  a  Gaussian,  5(r,  to)  is  reproduced  periodically  with 
a  period  Zp  -  n/a. 

It  would  appear  from  this  discussion  that  the  spectral  shift  should  follow  a 
periodic  evolution  pattern  with  period  Zp.  However,  this  is  not  the  case,  as  is  evident 
from  the  solid  curve  in  Fig.  6. 10  where  the  spectral  shift  is  plotted  as  a  function  of 
for  propagation  distances  covering  three  periods  (with  itoZp  ~  6600)  for  /totXi  =  20, 
=  10,  and  a  obtained  by  using  Eq.  (6.31).  Figure  6.10  is  drawn  for  a  fixed 
radial  distance  from  the  fiber  axis  (p  =  10)  rather  than  for  a  fixed  observation  angle. 
This  choice  is  made  because  for  large  propagation  distances  the  radial  distance  would 
exceed  the  fiber  dimensions  if  the  observation  angle  were  kept  fixed.  Figure  6.10  also 
shows  that  the  frequency  shift  indeed  becomes  zero  for  Zp  =  k/olq  [Oq  =  o:(coo)],  but 
its  maximum  and  ruinimum  values  become  larger  for  successive  periods. 
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Figure  6.10:  Frequency  shift  Av  as  a  function  of  propagation  distance  in  a  graded-index  fiber 
(solid  line).  The  frequency  shifts  are  calculated  for  observation  at  a  fixed  distance  IO/iIq  from  the  center 
of  the  fiber  and  ioOj  =  20  and  kQOg  =  10.  The  dashed  line  shows  the  frequency  shifts  when  the 
frequency  dependence  of  a  is  ignored  by  setting  [a(b)o)/llo  =  0.0(X)48]. 


The  physical  origin  of  the  non-periodic  nature  of  the  frequency  shift  can  be 
traced  back  to  the  dispersive  nature  of  the  graded-index  media  that  makes  a  frequency 
dependent  Indeed,  if  a  is  replaced  by  Oq  ,  we  obtain  the  behavior  indicated  by  the 
dashed  curve  in  Fig.  6.10.  It  is  clear  from  this  curve  that  when  the  frequency 
dependence  of  a  is  ignored,  the  frequency  shift  Av  shows  periodic  behavior  with 
period  Tt/o^  .  When  a  is  allowed  to  vary  with  frequency,  Av  becomes  non-periodic. 
This  feature  may  be  understood  by  noting  that  the  period  Zp  =  n/a  itself  becomes 
frequency  dependent  Since  the  argument  ocz  of  the  trigonometric  functions  appearing 
in  Eqs.  (6.17)  and  (6.19)  is  frequency-dependent  we  may  expect  z-dependent  changes 
in  the  frequency  shift  Thus  we  conclude  that  the  spectral  shifts  occurring  in  an 
inhomogeneous  medium  are  strongly  affected  by  the  dispersive  nature  of  the  mediunx 
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Appendix  C:  Derivation  of  the  expression  for  the  spectrum  of 
the  field  Eq.  (6.15) 

The  spectrum  of  the  field  at  any  point  r  is  given  by  the  expression 
S(r;Q))  =  ||A:*(r;p^;a))A'(r;p|2;ti))lV(rj,r2;(D)d^Pid^p(2  .  (Cl) 


Here  Ar(r;p;tD)  is  the  propagator,  given  by  Eqs.  (6.8),  and  is  the  cross- 

spectral  density  in  the  source  plane,  given  by  Eq.  (6.14).  On  substituting  from  Eq.  (6.14) 
into  Eq.  (C  1)  we  find  that 

5(r,co)  =  5<°>((0)f  --— -1  Bix,a>)B(y;Gy)  ,  (C  2) 

V27tsina2/ 


where 


Bix 


(C3) 


To  perform  the  two-dimensional  integration  we  introduce  the  average  and  difference 
variables 


Y2  =  ^2“^i- 
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(C  10) 


and  use  Eq.  (C  8)  again.  We  then  find 


ab 


ka 


:)■[- 


f  kacosoLz 

1^2flsinaz>/ r  [labimoLzJ 


(Cll) 


'Ve  note  that 


1- 


f  to  cos  02 
1 2aZ>sinaz  ) 


1 


(lab  sin  02) 


j|(2aZ>sinaz)^ -(itacostxz)'^j  ,  (C  12) 


which,  by  Eq.  (C  10)  can  be  simplified  to  the  fonn 


tocosaz  Y 
labsinazj 


(C13) 


On  substituting  from  Eq.  (C  11)  and  Eq.  (C  13)  into  Eq.  (C  2)  and  using  a  similar 
expression  for  B(y;  co)  we  obtain  the  expression 


5(r;to)  =  5(®^(©)^: 


)(  fexD. 

2  2 
X  +y 

ka 

r 

}\  1 
V2aZ>sinazy 
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4 

(C14) 


which  can  be  written  in  the  form 


5(r;a))  =  5(°>(tD) 


(C15) 
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A  =  2abCj 


sincxz 

02 


(C16) 
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Appendix  D:  Proof  of  the  periodic  reproduction  of  the 

source  spectrum 


When  the  propagation  distance  satisfies  the  condition  az  =  mn  for  a  positive 
integer  m,  some  of  the  factors  in  Eq.  (6.13)  become  singular.  However,  the  spectrum  of 
the  field  at  such  propagation  distances  remains  well  defined.  In  this  appendix  we  use  the 
method  of  stationary  phase  to  evaluate  the  spectrum  in  the  limit  as  z  ^  mn/a.  We  start 
with  Eq.  (6.13)  and  rewrite  it  as 


2 

S(r;tD)  =  jj  jj ^1  ^2 


(Dl) 


where 


sinoz  (D2) 

Since  A  — ♦  oo  as  z  mn/a,  we  can  evaluate  the  integral  using  the  method  of  stationary 
phase*.  According  to  this  method,  if 

/(A)  =  J/(r)exp[iA(l>(0]dr  ,  (D3) 


one  has,  under  fairly  general  conditions 

/(A)~cxp[iA(|)(d)]/(d)  C3^exp|^-^j  +  0(A"^^2)  as  A->oo.  (D4) 
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Here  d  is  a  zero  of  <|>'(0  [assuming  there  is  only  one]  and  p.  =  sgn[<|>"(d)].  /-pplying  these 
formulas  to  Eq.  (D  1),  with  the  definitions 

<(>(42)  =  “^^2]  ’  ^  5) 


and 


/(^2)  =  ^(^i’^2«'ni.'n2’“)  • 


(D6) 


we  find  that 


xexm 


f  x'^i 

i  cosaz  ' 

<2  cosozjj 

cosaz  ‘  \  f 


(D7) 


We  repeat  the  same  procedure  for  the  integration  over  the  variables  4i*  “Hi*  'n2 
asymptotic  limit  A  -> «».  All  the  phase  factors  cancel  and  the  final  result  is 

5(r;o))  =  wf-^,-^,-^.-^;o)l  .  (D8) 

V cosaz  cosaz  cosaz  cosaz  J 


Since  cosaz  =  x  1  we  see  from  the  last  equation  that  the  spectrum  of  the  source  is 
completely  reproduced  at  propagation  distances  z  =  Imnla.  For  propagation  distances 
z  =  (2m  +  l)7c/a  the  spectrum  of  the  source  is  again  reproduced  but  is  spatially  inverted. 
For  a  symmetric  intensity  profile  such  as  Gaussian,  the  source  spectrum  and  the  intensity 
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distribution  are  both  reproduced  at  z  =  mn/a.  It  can  be  shown  that  Eq.  (D  8)  is  exact  at 
z  =  mn/a  in  spite  of  our  use  of  an  asymptotic  approximation  for  evaluating  the  integrals  in 
Eq.  (D  1). 

1  N.  Bleistein  and  R.A.  Handelsman,  Asymptotic  Expansions  of  Integrals,  (Holt, 
Rinehart  and  Winston,  New  Yorit,  1975),  Chapter  6. 
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Summary 

In  this  thesis  we  investigated  effects  of  source  correlations  on  the  fields  radiated 
by  sources  of  any  state  of  coherence.  In  particular  we  studied  the  effects  of  spatial 
correlations  of  source  fluctuations  on  the  total  emitted  power  and  on  the  spatial  and 
spectral  distribution  and  the  spectrum  of  the  emitted  radiation. 

After  a  general  introduction  presented  in  Ch^ter  1  we  reviewed  some  elements 
of  coherence  theory  that  are  essential  to  the  discussion  in  this  thesis.  Chapter  2  consists 
of  two  parts.  In  the  first  part  we  introduce  the  concept  of  radiation  efficiency  and  we 
discuss  the  efficiency  of  planar,  secondary,  Gaussian  Schell-model  sources.  We  show 
that  for  sources  of  this  class  the  radiation  efficiency  increases  with  increasing  source 
size  and  with  increasing  spatial  correlation  length  of  the  source  fluctuations.  We  found 
that  when  the  linear  dimensions  of  the  source  exceed  several  wavelengths,  the  radiation 
efficiency  is  typically  over  90%,  irrespective  of  the  source  correlation  length^.  A 
different  result  is  obtained  for  the  radiation  efficiency  of  three-dimensional  primary 
sources;  in  this  case  we  find  that  the  radiation  efficiency  does  not  increase  with 
increasing  source  size  or  with  increasing  correlation  length.  In  fact,  we  find,  for 
example,  that  the  radiation  efficiency  of  a  uniform,  spherically  symmetric,  coherent  and 
co-phasal  source,  decreases  with  increasing  source  radius. 

For  a  primary  partially  coherent,  three-dimensional  sources  with  a  given 
intensity  profile  we  developed  a  method  which  makes  it  possible  to  determine  the 
spatial  correlation  function  that  maximizes  the  radiation  efficiency.  Using  this  method 
we  showed  that  the  radiation  efficiency  of  an  isotropic  quasi-homogeneous  source  is 
maximized  when  the  degree  of  spatial  coherence  is  given  by  smkr'Ikr'.  This  result  is 
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significant  in  that  it  shows  from  the  point  of  view  of  coherence  theory  that  the  spatial 
correlation  of  blackbody  radiators  gives  rise  to  the  maximum  radiation  efficiency  of 
isotropic  quasi-homogeneous  sources. 

The  discussion  of  the  radiation  efficiency  was  formulated  using  the  space- 
frequency  representation  and  the  results  therefore  apply  on  a  frequency  by  frequency 
basis.  In  chapter  3  we  focused  our  attention  on  the  frequency-dependence  of  fields 
generated  by  partially  coherent  sources  by  investigating  the  effects  of  source 
correlations  on  the  spectrum  of  the  radiation.  We  analyzed  a  basic  physical 
configuration  in  which  source  correlation  can  be  manifested,  i.e.,  that  of  two  small 
sources.  Using  this  simple  example  we  illustrated  possible  modifications  of  field 
spectra  depending  on  the  choice  of  the  correlation  function.  Our  calculations  included 
examples  of  line  narrowing,  line  broadening,  line  shifting  and  line  splitting. 

The  range  of  possible  modificaticm  of  field  spectra  can  be  significantly  enhanced 
by  using  a  system  which  contains  more  than  two  sources.  One  such  system  is  an  array 
of  2N  sources  which  we  also  analyzed  in  the  last  part  of  chapter  3. 

The  discussion  in  chapter  3  of  spectral  modulation  in  scalar  theory  was 
extended  in  chapter  4  to  the  full  electromagnetic  case.  We  considered  the  effect  of 
ccHTClation  on  the  spectrum  of  the  electromagnetic  field  produced  by  partially  correlated 
linear  dipoles.  In  this  chapter  we  were  mainly  concerned  with  the  dependence  of  the 
angular  distribution  of  the  radiant  intensity  and  with  the  shape  of  the  spectrum  in 
various  directions  of  observation.  We  showed  that  the  the  number  of  lobes  of  the 
radiant  intensity  produced  by  the  partially  correlated  dipoles  are  determined,  just  as  in 
the  deterministic  case,  mainly  by  the  separation  distance  between  the  two  dipoles, 
whereas  the  direction  and  the  shape  of  the  lobes  were  largely  determined  by  the 
correlation  between  the  two  dipoles. 
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In  the  analysis  of  chapters  3  and  4  we  assumed  that  we  can  choose  certain 
forms  of  spatial  correladons.  In  chapter  5  we  analyzed  the  radiation  from  a  source  that 
consisted  of  a  finite  number  of  coherent  modes.  The  spatial  correlation  of  this  soun  e 
then  depends  both  on  the  number  of  the  modes  and  on  their  relative  strengths.  For 
sources  that  consist  of  Hermite-Gaussian  modes  we  showed  the  spectral  effects  in  the 
far  and  in  the  near  zone  and  dieir  dependence  on  the  mode  composition  of  the  source. 
We  also  examined,  for  comparison,  the  small  spectral  shifts  that  occur  on  propagation 
from  sources  that  consist  of  a  single  coherent  mode.  In  the  analysis  of  chapter  5  we 
introduced  the  concept  of  "spectral  modifier"  which  is  a  correlation-dependent  factor 
that  determines  the  spectral  changes  in  the  emitted  radiation. 

Chapters  2-5  pertain  to  propagation  in  free  spac''.  In  chapter  6  we  considered 
propagation  of  partially  coherent  tight  in  homogeneous  and  inhonx)geneous  media.  We 
also  studied  how  the  spectrum  changes  as  a  function  of  the  propagation  distance  from 
the  source  plane.  Our  results  show  that  spectral  shifts  are  en  ced  in  a  medium  whose 
index  of  refraction  is  larger  than  unit^.  Within  a  given  homogeneous  medium  tl.«’ 
spectral  shixes  usually  develop  rapidly  with  increasing  propagation  distance  from  the 
source  plane.  In  the  case  of  a  non-dispersive  graded-index  fiber,  we  showed  that  the 
spectral  shifts  are  periodic  along  the  propagation  diicction.  When  dispersion  effects  are 
taken  into  account  we  showed  that  there  exists  deviation  of  the  spectral  shifts  f.om  the 
periodic  form. 

The  research  described  in  tliis  thesis  '  conducted  between  1986  and  1989. 
During  this  period  many  other  contributions  were  made  in  t*'is  field  by  numerous 
authors.  It  seei.is  that  there  is  a  great  potential  for  further  work  in  this  area,  especially 
in  connection  with  astrophysics,  spectruradiometry  and  scattering. 
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*  This  statement  assumes  that  the  correlation  length  is  not  less  than  that  of  a 
thermal  source. 
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